AD-A268  521 


#• 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey ,  CaMenia 


THESIS 


OTlC 
I  E  ^ 


NUMERICALLY  SOLVING  A  TRANSIENT  HEAT  CONDUCTION 
PROBLEM  WITH  CONVECTION  AND  RADIATION 


David  J.  Albert 


June  1993 


Thesis  Advisor 


Jeffery  Leader 


Approved  for  public  release;  distribution  Is  unlimited. 


3  8  24  05.5. 


93-19719 

lillllil 


REPORT  DOCUMENTATION  PAGE 


form  Approvtrt 
0MB  No  0704  01BB 


»UW.C  >nionir<9  l>u>d«n  <o»  iht»  eoll«tw  of  iiHwirwiion  .«  iq  «.««!«  i  how  on  >«tpo<nc.  fn(lud.nq  (h«  imt  lot  •fvtirwina  ••niructiom  emlino  data  u>urc« 

<0«M1<W.  of  inform.tio«,  .nclud<i>9  wffniw  fm  rMuonf  tM*  Ourdw  to  Mmhinqtoa  Htooouon,,,  t«r..,«,.  OtrKloroi*  Oo*f«itom  onO  MmLIIltiVHMrM 

0«..t  s»..t«  U0«  a,i.«9I0«.  va  HiOiJSOi  ,nd  to  ttw  Offtc*  of  i^onogtmtot  to*  •udq*,.  »w,»o.ii  Hoduolon  rtofKI  oTjOSO J 


1.  AGINCV  USE  ONLY  (Lttvt  bistyk) 


2.  REPOKT  DATE 

June  1993 


i.  REfOAT  TYRE  AND  OATES  COVERED 
Master's  Thesis 


A.  TITLE  AND  SUSTITLE 

MUMERICALLY  SOLVING  A  TRANSIENT  BEAT  CONDUCTION  RROBLEN  NITH 
CONVECTION  AND  RADIATION 


«  AUTHOR(Sj  — — 

David  J.  Albert 


7.  PERFORMING  ORGANIZATION  NAME(S)  ANO  AOORESS(ES) 
Naval  Postgraduate  School 
Monterey,  CA  93943-5000 


PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  ANO  AOORESS<ES) 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 


11.  SUPPLEMENTARY  NOTES  — - - 

The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect 

the  official  policy  or  position  of  the  Department  of  Defense  or  the  U.S.  Government, 


IJl.  DISTRIBUTION /AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited 


IZb.  DISTRIBUTION  CODE 
A 


13.  AtSTRACT  (Afjjiimum  200  wordi) 

.  .  j ^5*”®!.^.*”*^  aurfacB  temporaturs  dlatrlbution  is 

*p*'®*^*  subjected  to  coollng 
^  convection  and  tadlatlon.  in  the  study,  the 

V  '3?“*'''  Vl"®  problem  Is  reduced  to  a  singula? 
nonlinear  Volterra  Integral  equation  of  the  second  kind  Lino 
the  Integral  transform  method.  Several  numerical  technlaues 

of%he*orobfem.  ^  ®"  approximate  solution 

RunqS%SSa^  meVh ®  esslve  approximations,  the 

"fthod,  and  the  finite  difference  method.  The 
equation  is  solved  numerically  by  the  Runge-Kutta 

--i®®  problem,  and  the  solutions  are  compared  with  those 
generated  by  the  Runge-Kutta  method.  All  the  numerical 

difficulties®  Limitations  and 

difficulties  Involved  In  these  echemes  are  discussed.  At  the 

??opos?d."'*"**^^‘="'  algorithm  for  solving  the  problJi  1? 

U.  SUBJECT  TERMS  '  '  "  '  "  '  ~  I  fl  flUlurff I  TF 

Numerical  Analysis,  Heat  Equation,  Runge-Kutta,  Finite 

Difference,  Volterra  Integral  Equation  ■■  .-IlOQ— e— 

IB.  PRICE  CODE 


17.  security  ClASSeiCATION 

11  SECURITY  CLASSIFICATION 

19.  SECURITY  CLASSIFICATibN 

OF  REPORT 

OF  THIS  PAGE 

OF  ABSTRACT 

Unclass 

Unclass 

Unclass 

Approved  for  public  release;  distribution  is  unlimited 


Numerically  Solving  A  Transient  Heat  Conduction  Problem 

With 

Convection  and  Radiation 


by 

David  J.  Albert 

Lieutenant  ,  United  States  Navy 
B.S.,  University  of  North  Carolina,  1985 

Submitted  in  partial  fulfillment 
of  the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  APPLIED  MATHEMATICS 

from  the 


NAVAL  POSTGRADUATE  SCHOOL 


ABSTRACT 


The  transient  surface  temperature  distribution  is 
determined  for  the  flat  plate  and  sphere  subjected  to  cooling 
by  combined  convection  and  radiation.  In  the  study,  the 
initial  boundary  value  problem  is  reduced  to  a  singular 
nonlinear  Volterra  integral  equation  of  the  second  kind  using 
the  integral  transform  method.  Several  numerical  technic[ues 
are  introduced  in  an  attempt  to  find  an  approximate  solution 
of  the  problem:  The  method  of  successive  approximations,  the 
Runge-Kutta  method,  and  the  finite  difference  method.  The 
integral  equation  is  solved  niimerically  by  the  Runge-Kutta 
method  of  orders  1,  3,  and  5.  In  addition,  the  finite 
difference  method  is  implemented  to  solve  the  initial  boundary 
value  problem,  and  the  solutions  are  compared  with  those 
generated  by  the  Runge-Kutta  method.  All  the  nximerical 
results  are  presented  graphically.  Limitations  and 
difficulties  involved  in  these  schemes  are  discussed.  At  the 


end,  a  ntimerical  algorithm  for  solving  the  problem  is 
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I.  ANALYTIC  SOLUTIONS  OF  THE  HEAT  EQUATION  SUBJECT  TO 
CONVECTIVE  AND  RADIATIVE  BOUNDARY  CONDITIONS 

A.  INTRODUCTION 

During  the  60 's,  space  technology  advanced  so  much  that 
the  research  of  the  temperature  behaviour  of  bodies  exposed  to 
a  deep  space  environment  became  crucial.  In  particular, 
transient  heat  or  cooling  of  solids  of  different  shapes  by 
convection  and  thermal  radiation  was  becoming  highly  important 
in  many  engineering  applications.  An  example  of  these 
applications  is  the  temperature  distributions  of  rocket 
motors.  An  extensive  investigation  of  the  problem  has  been 
conducted  and  a  lot  of  literature  on  the  subject  was  published 
during  the  60 's  and  70 's.  A  detailed  review  of  most  of  these 
papers  Is  not  intended  here;  instead  a  brief  summary  of  the 
major  ones  will  be  given. 

As  early  as  1962,  Fairall,  et.al.[6]  generated  a  numerical 
solution  for  the  problem  using  an  explicit  finite  difference 
scheme;  this  paper  served  as  pioneer  work  in  the  area  of  the 
research.  Later,  various  finite  difference  schemes  were 
devised  to  deal  with  the  nonlinear  boundary  condition.  The 
main  difficulty  in  these  schemes  is  the  appearance  of  severe 
oscillations  in  the  determined  temperature  values  for  high 
heat  flux  situations.  Von  Rosenberg  [10]  proposed  a  hybrid  of 
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an  iterative  technique  and  implicit  finite  difference  schemes 
to  deal  with  the  nonlinear  boundary  condition.  On  the  other 
hand,  Crosbie  and  Viskanta  [3,4]  transformed  the  governing 
equations  into  a  nonlinear  Vol terra  integral  equation  of  the 
second  kind  and  applied  the  method  of  successive 
approximations  to  solve  the  integral  equation.  Milton  and 
Goss  [8,9]  developed  some  heuristic  stability  criteria  for 
explicit  finite  difference  schemes  with  nonlinear  boundary 
conditions.  It  turns  out  that  a  very  restrictive  time  step  is 
required  for  numerical  stability  which  may  result  in  requiring 
a  prohibitive  amount  of  computer  time  to  calculate  the  long 
time  evolution  of  the  solutions.  Williams  and  Curry  [12] 
surveyed  several  methods  for  treating  the  nonlinear  boundary 
condition  in  implicit  schemes  and  compared  their  accuracy  and 
efficiency. 

Nonlinearity  is  commonplace  in  natural  phenomena. 
Unfortunately,  a  nonlinear  problem  often  doesn’t  lend  itself 
to  a  closed  form  solution.  The  problem  of  transient  heat- 
conduction  in  a  solid  becomes  nonlinear  when  the  surface  of 
the  body  is  subjected  to  thermal  radiation.  When  energy 
transfers  through  the  wall  of  a  body,  two  cases  arise: 
convection  and  thermal  radiation.  The  convective  heat 
transfer  describes  the  situation  where  heat  is  dissipated 
according  to  Newton's  Law  of  Cooling,  which  states  that  the 
rate  at  which  heat  is  transferred  from  the  body  to  a 
surrounding  is  proportional  to  the  difference  in  temperature 
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between  the  body  and  the  environment.  The  boundary  condition 
that  describes  convection  is  nonlinear  except  for  the  case 
where  the  heat-transfer  coefficient  is  independent  of  surface 
temperature,  which  is  technically  called  forced  convection. 
The  radiative  heat  transfer  is  based  on  the  Stefan-Boltzmann 
Law,  which  states  that  the  heat  flux  is  proportional  to  the 
difference  between  the  surface  temperature  to  the  fourth 
power  and  the  source  temperature.  Pure  radiation  or  pure 
convection  occur  whenever  one  mode  of  energy  transfer 
predominates  over  the  other. 

It  is  the  purpose  of  this  thesis  to  consider  the  one¬ 
dimensional  transient  heat  conduction  problem  resulting  from 
a  combined  convective  and  radiative  heat  flux  with  the 
objective  of  determining  the  surface  temperature  fields  using 
the  numerical  methods  which  are  discussed  in  this  study. 
Another  purpose  of  this  thesis  is  to  explore  the  limitations 
and  difficulties  involved  in  these  schemes.  References  to  the 
work  done  in  similar  areas  are  presented  to  allow  the  reader 
further  investigation. 

Analytic  solutions  are  derived  in  one  dimension.  However, 
the  resulting  solutions  are  not  in  closed  form,  and  thus 
impractical  to  use.  Hence,  numerical  techniques  will  be 
studied  and  employed  in  the  computer  in  an  attempt  to  find  an 
approximate  solution.  Numerical  results,  found  by 
implementing  some  of  the  numerical  methods  discussed  below, 
will  be  presented  and  compared.  In  the  conclusion,  a 
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numerical  scheme  is  proposed  as  an  alternative  to  the  existing 
methods.  It  is  open  to  the  readers  for  justification. 

Sections  1(C)  and  1(D)  describe  the  derivation  of  the 
integral  representations  of  the  one  dimensional  transient  heat 
conduction  problem  subjected  to  a  combined  convective  and 
radiative  boundary  condition  in  a  rectangular  coordinate 
system.  Two  integral  transform  methods,  namely  the  Laplace 
transform  and  the  eigenvalue  exp-ansion,  are  presented. 
Observation  and  comparison  are  made  for  the  integral  equations 
to  yield  some  useful  information  about  the  solutions. 

In  Chapters  II  and  III,  numerical  methods  for  the 
solutions  of  the  nonlinear  Volterra  integral  equations  of  the 
second  kind  are  described.  In  particular,  the  method  of 
successive  approximations  and  the  Runge-Kutta  method  are 
outlined  in  detail.  A  brief  remark  is  given  for  their 
advantages  and  limitations  in  finding  solutions  to  the 
integral  equation. 

Chapter  IV  describes  a  numerical  method  which  is  directly 
applied  to  the  governing  partial  differential  equation.  The 
technique  is  called  the  finite  difference  method.  It  is 
basically  a  hybrid  of  finite  difference  techniques  and  an 
iterative  scheme  proposed.  A  suggestion  is  made  for  the 
improvement  of  the  algorithm. 

In  Chapter  V,  numerical  results  produced  by  some  of  the 
discussed  numerical  schemes  are  presented.  The  implementation 
of  various  methods  gave  a  practical  sense  of  their  advantages 
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and  limitations.  Graphs  and  tables  are  set  up  in  such  a  way 
that  a  comparison  can  be  made. 

In  the  next  section,  a  statement  of  the  problem  is  given. 
In  the  statement,  the  basic  assumptions,  the  governing 
equation  and  the  boundary-initial  conditions  are  included. 

B.  STATEMENT  OF  THE  PROBLEM  FOR  OBTAZMIMO  THE  SURFACE 
TEMPERATURE 

Considering  the  one-dimensional,  transient,  conduction 
heat  transfer  problem  with  combined  convection  and  radiation 
at  its  surface,  the  following  assumptions  have  been  made: 

1.  One-dimensional  heat  transfer  to  a  solid  of  a  finite 
length. 

2.  The  solid  medium  is  pure,  isotropic,  homogeneous,  and 
opaque  to  thermal  radiation. 

3.  All  thermodynamic  and  transport  properties  are 
independent  of  temperature. 

4.  The  solid  does  not  contain  any  heat  sources  or  sinks. 

5.  The  fluid  is  transparent  to  thermal  radiation. 

6.  The  fluid  temperature  and  the  ambient  temperature  are 
constant . 

The  non-dimensional  form  of  the  governing  partial 
differential  equation  for  the  temperature  U(x,t)  and  the 
appropriate  initial  boundary  conditions  are 


d^U  dU 


0<X<1, 


t>0; 


(1.1) 
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with  initial  condition 


U(x,0)  =  g(x) 


(1.2a) 


and  boundaxy  conditions 

(i)  -  ciiUiO.t)  =  0  (1.2jb) 

(ii)  -  03C;(l,t)  =  -hC7Ul,t)  .  (1.2C) 

Note:  ffi  and  CI2  can  be  any  real  number,  except  both  cannot  be 
zero  at  the  same  time.  03  is  a  non-zero  real  number,  and  h  is 
a  positive  real  number. 

The  next  section  will  deal  with  solving  the  partial 
differential  equation  (1>1)  with  initial  and  boundary 
conditions  (1.2a-c)  by  the  Laplace  transform  method  and  the 
eigenvalue  expansion  method.  As  an  illustration,  two  special 
cases  with  specific  values  of  Cj,  ttj,  a3,  and  h  will  be 
considered,  and  the  analytic  solutions  of  these  cases  at  the 
surface  will  be  derived.  It  will  be  shown  that  the  surface 
temperature  satisfies  a  singular  Volterra  integral  equation  of 
the  second  kind.  At  the  end  of  the  chapter,  we  will  present 
the  solutions  and  indicate  some  useful  information  about  the 
integral  equations. 
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C.  THB  LAPLACE  TRAM8FOSM  METHOD 


In  this  section,  the  Laplace  transform  of  equation  (1.1) 
with  associated  boundary  conditions  (1.2b,c}  is  first  obtained 
with  respect  to  time.  The  resulting  boundary  value  problem  is 
in  terms  of  the  Laplace  transform  of  the  required  solution. 
Next,  the  ecpiations  are  solved  for  the  transformed 
temperature,  and  the  solution  of  the  stated  problem  can  be 
found  by  taking  the  Inverse  Laplace  transform  of  the 
transformed  solution.  From  experience,  it  can  be  expected, 
the  Laplace  inversion  is  of  some  difficulty.  To  simplify  the 
situation,  specific  values  of  Sx,  Sj,  a^,  and  h  are  considered 
so  that  the  inverse  process  is  practical  without  loss  of 
generality.  It  should  be  noted  that  there  does  exist  an 
inverse  Laplace  transform  for  other  cases  of  a  more  general 
nature . 

Now,  define  the  transform  of  the  temperature  function, 
U(x,t),  with  respect  to  time  as  follows 


aE[  C7(x,  t)  ]  (s)  =  J  u(x,  t)  e'*'dt  =  C7(x,s). 


(1.3) 


After  the  transformation,  the  temperature  function  becomes  a 
function  not  only  of  x  but  also  of  the  parameter  s.  Assuming 
that  the  derivatives  with  respect  to  x  pass  through  the 
transform  (differentiation  can  be  accomplished  before 
integration) ,  we  have 
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(1.4) 


y. au(x,t),,.,  .  p3u(x,t)^...^^ .  ag(x,s) 
OX  /O  OX  Ox 


gj[  (^)  ^  pyu(x,  t)  e-stgffc  =  (1.5) 

Ox^  /o  ox^  Ox^ 

The  rule  for  transforming  a  derivative  with  respect  to  time 
can  be  found  using  integration  by  parts.  Thus,  the  Laplace 
transform  of  the  derivatives  of  U(x,t)  with  respect  to  the 
transformed  variable  t  is  given  by 

]  (s)  =  =  sU(x,s)  -  U{x,0) .  (1.6) 

Now,  applying  the  Laplace  transform  to  the  initial -boundary 
value  problem  (1.1),  (1.2a-c)  we  remove  all  time  derivatives. 
Holding  s  fixed,  we  have  the  following  ordinary  differential 
equation  in  x 


0<x<l 


(1.7) 


with  boundary  conditions 

tt^— -  a2l/(0,s)  «  0,  for  X  =  0  (1.8a) 
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~  *  -hgt[U*(l,s)],  for  x^  1  (1.8b) 

Notice  that  the  initial  condition,  g(x),  is  incorporated  in 
the  ordinary  differential  equation.  In  order  to  solve  (1.7) 
and  (1.8a,b),  we  must  first  solve  for  the  general  solution  of 
the  corresponding  homogeneous  differential  equation  and  a 
particular  solution  of  (1.7)  satisfying  (1.8a,b).  Now, 
consider  the  general  solution  of  the  homogeneous  equation  for 
(1.7), 


UHcm(x,s)  =  +  Be^% 


(1.9) 


where 


^1,2  = 


(1.10) 


which  are  given  by  the  roots  of  the  auxiliary  equation 


-  s  =  0. 


(1.11) 


In  the  following  paragraph  we  employ  the  method  of  variation 
of  parameters  to  solve  for  a  particular  solution  of  (1.7). 
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Let 


Up(x,s)  =  C/iVi(x,s)  +C4v2(x,s),  (1.12) 

be  a  paxrticular  solution  where  Ui(x,s)  and  U2(x,s)  are  any  two 
linearly  independent  solutions  of  the  corresponding 
homogeneous  equation.  In  this  case,  choose  Ui(x,s)  »  e^“  and 
U2(x,s)  *  The  object  here  is  to  find  Vi(x,s)  and  V2(x,s) 

such  that  the  following  equations  are  satisfied 

eV«vi{x,s)  +  e-'^v^(x,s)  =0,  (1.13) 

(X, s)  -  ^e‘^^*V2(x,s)  =-  g(x)  (1.14) 


By  Cramer's  rule, 


v'i(x,  s)  = 


g{x)  e-^ 
-2y/S 


(1.15) 


and 


(x,s)  . 

2yfS 


(1.16) 


By  integrating  (1.15)  and  (1.16),  we  obtain 


Jo  2^ 


(1.17) 
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V,U,S)  .  -v,(0,s) 

Jo  2^ 


(1.18) 


Thus,  the  general  solution  to  (1«7)  and  (1.8a,b)  is 

U(X,S)  =  U^{X,3)  +  Uj,{X,8)  ll.iyj 

that  is, 

U{x,s)  *  Ae^  +  Be~^  +  e'^VjCx,  s)  +  e‘'^V2(x,  s)  ,  (1.20) 

where  A  and  B  are  arbitrary  constants  and  Ui(x,s),  U2(x,s)  are 
given  by  (1.17)  and  (1.18),  respectively.  To  determine  A  and 
B,  boundaxry  conditions  (1.8a,b)  are  used  along  with  the 
following  procedure.  The  derivative  of  U(x,s)  from  equation 
(1.20)  is  found  to  be 

=  A/se'^  -  Byfse''^  +  (x,  s)  + 

e''**Vi(x,s)  +  e"'^V2{x, s)  -  (x,  s)  .  (1.21) 

Let  X  =  0.  Then  (1.15)  and  (1.16)  give 

Vi(0,S)  +  V2(0,S)  =  0  . 

(1.18a),  (1.20),  and  (1.21)  then  imply 
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Oi  [AV^  -  s)  -^/SVjCO,^)] 

-  ttjEA  +  B  +Vi(0,s)  +V2(0,s)]  =0  (1.22) 

By  rearranging  the  tezrms,  (1.22)  becomes 
A(ai\/s  -  ttj)  -  B(ai\/s  +  ttj)  * 

+  a2)v2(0,s)  -  -  a2)Vi(0,s)  .  (1.23) 

similarly  let  x  =  1.  Then  (1.15)  and  (1.16)  give 

e'^l(l,s)  +  e"'^2(l/S)  *0  .  (1.24) 

Therefore  (1.8a),  (1.20),  and  (1.21)  then  imply 

Ayfae^  -  B^e'^  +  y/se'^^{X,s)  -  " 

ttjEAe'^  +  Be~^  +  e'^i(l,s)  +  e''^V2(l,s)]  = 

-  hS£[cr«(l,  t)]  .  (1.25) 

By  a  similar  manipulation  of  the  terms,  (1.25)  becomes 
A{yfse^  -  ttj©'^)  -  B{^/se~'^  +  = 

(v'se"'^  +  a3©‘'^)v2(l,s)  -  i^e'^  -  Sj©'^) (i, s)  - 

ASf[£7Ml,  t)3  .  (1.26) 
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Equations  (1.23)  and  (1.26)  form  a  system  of  two  equations  in 
the  two  unknowns  A  and  B.  By  Cramer's  rule,  A  and  B  are  as 
follows 


where 


,  _  numi 

-ds- 


(1.27) 


numl  - 

Vi  (0,  s)  +  (ttiv's+aj)  (Vj  (1,  s)  -Vj  (0,  s)  ]  }(^e'v^+a3e''^) 

-  [(aiVs+a2)Vi{l,s)]  (l,  t)  ]  (Oi/i+aj)  , 


and 


den  = 


(ttj  ’  a3ai)v^(e'^  +  e"'^)  +  ia^s  -  a^ctz)  [e^-e'^)  . 


_  num2 
den 


(1.28) 


where 


nim2  = 

((ttiVs-aj)  [vi(0,s)  -Vi(l,s)]  -(aj/s+a2)v2(0,s)}(v^e'^-a3e'^) 
-  Ii^[U*  (1,  t)]  {a^y^s-az)  *  [  (aiVs-a2)  Vjd,^)  ]  {^e'^^^a^e''^) 


and  den  is  the  same. 
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Thus,  the  general  solution  of  (1.7),  (1.8a,b)  is  given  by 
(1.20)  where  A,  B,  Vi(x,s),  and  U2(x,s)  are  given  by  (1.27), 
(1.28),  (1.17),  and  (1.18),  respectively.  Theoretically,  the 
analytic  solution  of  partial  differential  equation  (1.1)  with 
initial  and  boundary  conditions  (1.2a,b,c)  can  be  obtained  by 
taking  the  Laplace  inversion  of  U(x,s),  and  thus,  the  surface 
solution  can  be  found  by  putting  x  equal  to  one  in  U(x,s) .  In 
practice,  however,  the  inverse  Laplace  transformation  process 
is  highly  unstable  in  that  singularities  may  exist.  Also,  the 
transforms  are  difficult  to  find.  In  the  next  paragraph  below 
we  consider  two  special  cases  where  the  inversion  is  feasible. 
In  each  case,  values  for  the  parameters  correspond  to  a 
specific  geometrical  configuration  of  a  body. 

Casg  li  0i  =  1,  ttz  *  0,  aj  =  -1,  h  =  1 

This  set  of  values  corresponds  to  a  "flat  plate"  with  a 
given  initial  temperature  and  which  is  being  heated  or  cooled 
by  combined  convection  and  radiation.  The  term  "flat  plate" 
is  taken  here  to  mean  a  solid  slab  of  finite  thickness  which 
is  bounded  by  a  pair  of  vertical  lines  at  ±  ^  thus  having  a 
width  of  1.  Substituting  the  given  values  for  the  parameters 
in  (1.20)  we  obtain  the  transformed  surface  temperature 

U{l,s)  -  Ae'^  +  Be"^  +  s)  +  (l,  s)  ,  (1.29) 

where , 
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(1.30) 


[%^i(0,s)  +v^a(l,s)  -  %^2(0,s)  1  (^e~^  -  e'>^) 
/s(e'^  +  e'^)  +  s{e^  -  e~'^) 

^  -  [v/gVi(l,g)]  (y^ev^  +  e>^)  -  A£t[ty*{l/ 1:)]%/5 
v/s(e^  +  e"^)  +  s(e'^  -  e"^) 


j  ^  [^i(0,s)  -v^i(l,s)  -  Vsva(0,s)]  (^/se^  ^•  ev^) 

(aj  -  a3ai)v^(e'^  +  e''^)  +  (a^s  -  ajttj)  (ev^  -  e"'^) 

^  -  tiSSi[U*{l,t)\^/s  *  [^2(l,s)l  (v^e~^  -  e~v^) 

(02  -  03ai)v'5(e'^  +  e"^)  +  {a^s  -  a^cL^)  {e^  -  e'v^) 


,,(i,s)  . -f‘j£Me2±.'dx'-v.(0,s) 
Jo 

and 

VjCl,^)  =  p.2lifE_L£2C_obc'  +  V2(0,s) 
Jo  2^ 


Now  substituting  (1. 30) -(1. 33)  into  (1.29)  and  simplifying  the 
results  gives: 

rr/.  X  _  /5sr(x')e^'dx' + /5g(x0e-^W 
XJ  \X  g  S)  "•  '  '  ■"  — —  . .  -  —  N 

(e>/s  +  e"'^)  +  -  e~'^) 


+ 


-  tM[U*{l,  t)l  {e^  e~v^) 

(ev^  +  e'^)  +  v's(e'^  -  e~'^) 


(1.34) 
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Suppose  the  initial  temperature  is  1,  that  is, 

g(x)  -  1. 


(1.35) 


The  boundary  conditions  associated  with  the  given  values  of 
the  parameters  and  initial  condition  (1.35)  constitute  a 
cooling  process.  With  (1.35),  the  transformed  surface 
temperature  becomes 

_L  (eV2  _  e-Vi)  -  h^[U*{l,t)]  *  e-^) 

U{l,s)  =  - - - - - - -  .  (1.36) 

(e\/3  +  e"'^)  +  yfs{ey/s  -  e"'^) 

If  (1.36)  is  multiplied  through  by 

1  ^  *  e'^  ^(e^  -  e'^)  ^ 

Vs  (e'^  -  e“'^) 


and  then  simplified. 


U{l,s) 


1  _  yt[hU*{l,t)  +  U(l,  t)]  {e^  +  e-^) 
®  VsKe'^  -  e'^) 


(1.37) 


is  obtained.  Equation  (1.37)  is  ready  to  be  inverted.  In 
order  to  perform  the  inversion  of  (1.37),  the  following  two 
Laplace  transforms  have  to  be  computed 


S£-M-]  and 

3 


C£-l[ 


(gy/i  ^  e~^)  j 
^(gVi  -  e"'^) 
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In  fact  the  transforms  can  be  found  from  any  standard  Laplace 
transform  table.  By  the  convolution  theorem,  the  surface 
temperature  in  time  t  is  given  by 


U{1,  t)  =  1 


Jo 


[1 


Jr-1 


e 


^1t[C  -  X) 


+  :7(l,T)]dt  (1.38: 


and  by  the  Poisson  summation  formula  [14],  (1.38)  can  be 

written  as 

Uil.O  =  1  -  [t;Ml.t)+U{l,T)]  eft  .  (1.39) 


Hence,  the  problem  of  transient  cooling  of  a  flat  plate  by 
combined  convection  and  thermal  radiation  has  been  reduced  to 
solving  a  nonlinear  Volterra  integral  equation  of  the  second 
kind. 

Case  2;  tti  =  0,  Oj  =  -1,  Sj  *  1,  h  =  1 

This  set  of  values  corresponds  to  the  case  where  a 
spherical  body  of  radius  1  with  a  given  initial  temperature 
is  being  heated  or  cooled  by  combined  convection  and 
radiation.  Since  the  procedures  used  to  solve  the  problem  are 
basically  those  described  in  case  1,  the  mathematical  details 
will  be  omitted  and  only  the  main  steps  will  be  presented. 
Consider  equation  (1.20) ,  the  general  solution  of  the  boundary 
value  problem.  The  given  values  for  the  parameters  are  first 
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substituted  Into  (1.17),  (1.18),  (1.27),  and  (1.28).  Then, 

(1.20)  is  simplified  as  in  the  previous  case.  After  a  tedious 
calculation,  the  transformed  surface  temperature  is  given  by 

f^Sr(x')  e-^'dx'  -  f^ff(x')  e^'dx'  + 

U{l,s)  =  ^ ^ - 

(ev^  -  e~^)  -  +  e*'^) 

^  h  ^[U*{l,t)]  ie^  +  e-^) 

(ev^  -  e-'^)  -  /s{e>^  +  e'^)  ' 

Suppose  the  initial  temperature  is  chosen  to  be 

g(x)  =  X. 

Boundary  conditions  associated  with  the  given  values  of  the 
parameters  and  initial  condition  (1.41)  again  constitute  a 
cooling  process.  With  (1.41),  the  transformed  surface 
temperature  becomes 

17(1,  s)  =  —  -  — A  S£[t7*(l,  t)  ]  (e^  -  e  —  (1.42) 


(1.40) 


(1.41) 


which  is  now  ready  to  be  inverted.  In  order  to  perform  the 
Laplace  inversion  of  (1.42)  the  following  two  inverse  Laplace 
transforms  need  to  be  computed 


— ]  and 
s 


ae-M - rr 


(e^  - 
+  e"^) 


e-v^) 

-  (e'^ 


e-^) 


]  . 
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The  first  inverse  Laplace  transform  is  obvious.  However,  the 
second  one  is  not  so  obvious.  Details  of  the  derivation  of 
the  second  inverse  Laplace  transform  are  given  in  [ 1 ] .  The 
surface  temperature  in  time  t,  obtained  by  inverting  (1.42), 
is 

U{l,t)  =1  (1.43) 

where  is  k—  positive  root  of  the  transcendental 

equation 

Pk  *  tan  Pk  .  (1.44) 

Hence,  the  problem  of  transient  cooling  of  a  sphere  by 
combined  convection  and  thermal  radiation  has  been  again 
reduced  to  solving  a  nonlinear  Volterra  integral  equation  of 
the  second  kind.  As  we  have  mentioned  above,  one  of  the 
drawbacks  of  the  Laplace  transform  method  is  that  there  are 
only  a  few  cases  in  which  the  transformed  solution  can  be 
practically  inverted  into  the  required  solution.  In  the  next 
section,  the  eigenvalue  expansion  method  is  introduced  as  an 
alternative  to  the  above  method.  One  may  find  the  eigenvalue 
method  more  practical  for  solving  for  the  analytic  solution  of 
the  heat  equation  with  nonlinear  boundary  conditions. 

D.  THE  EIGENVALUE  EXPANSION  METHOD 

The  fundamental  idea  of  the  eigenvalue  expansion  method  is 
to  transform  the  given  boundary  value  problem  by  the 
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eigenfunctions  obtained  from  the  associated  eigenfunction 
problem.  By  the  completeness  theorem  (which  states  that  any 
piecewise  smooth  function  can  be  represented  by  a  generalized 
series  of  eigenfunctions)  we  can  show  that  separation  of 
variables,  i.e.,  u(x,t)  =  X(x)T(t) ,  may  lead  to  the  solution 
of  the  problem  expressed  as  an  infinite  sum  of  the 
eigenfunctions  with  appropriate  coefficients  determined  by  the 
orthogonality  property  of  eigenfunctions.  Applying  these 
procedures  to  the  partial  differential  equation  (1.1)  and 
initial  boundary  conditions  (1.2a,b,c)  yields  the  following 
main  results 

+  p2jf(x)  =0,  0  <  X  <  1  (1.45) 

dx^ 


with  boundary  conditions 


-  a^XiO)  =  0, 


and 


-  ajJfd)  =  0 
dx  ■* 


(1.46) 


(1.47) 


Parameters  ai  and  a2  can  be  any  real  number  except  they  cannot 
be  zero  at  the  same  time.  is  a  non-zero  real  number. 
According  to  the  theory  of  ordinary  differential  equations, 
the  general  solution  of  (1.45)  is 
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X(x)  =  CiCOs(Pjc)  +  CjSinCPx)  . 


(1.48) 


Applying  boundary  conditions  (1.46)  and  (1.47)  to  equation 
(1.48)  gives  the  following  system  of  equations 

ttipCa  =  ttjCi  (1.49) 

(CjP  -  Cia3)cosP  =  (CiP  +  C2a3)sinp  .  (1.50) 

Note  that  boundary  value  problem  (1.45  -  1.47)  is  in  the 
class  of  Sturm-Liouville  problems  for  which  all  eigenvalues 
are  real  and  the  eigenfunctions  corresponding  to  different 
eigenvalues  are  orthogonal.  Thus,  if  the  parameters  in  (1.49) 
and  (1.50)  are  specified,  there  will  exist  eigenvalues,  p„, 
where  n  »  1,2,...,  and  the  corresponding  eigenfunctions, 
Xn(x),  such  that  the  temperature  function,  U(x,t),  can  be 
expanded  in  a  Fourier  expansion  of  the  form 


U{x,t)=YZ^^Ujt)Xjx),  (1.51) 


where  the  Fourier  coefficients,  U„(t),  are  given  by 

U„{t)  =  f^U{x,C)  X„{x)dx  .  (1.52) 
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Now,  taking  the  finite  Fourier  integral  transform  of  the  heat 
equation  (1.1)  with  respect  to  X„(x)  gives 


t)  Ar„(x)  dx  =  /„‘0A-„(x)  dx  . 


(1.53) 


Performing  integration  by  parts  of  the  right  hand  expression 
in  equation  (1.53)  and  substituting  in  (1.52)  yields  the 
following  ordinary  differential  equation  for  U„(t) 


dUjC) 

dt 


dU{l,  t) 


^„(1) 


-U(l,t)xUl) 


+  U{0,t)xUo)  +  ru(x.t)x"(x)  dx  . 

«^0 


(1.54) 


With  boundary  conditions  (1.46)  and  (1.47),  the  right  hand 
side  of  (1.54)  can  be  simplified.  Then,  by  the  integrating 
factor  method,  the  solution  of  equation  (1.54)  can  be 
obtained.  Hence,  the  resulting  integral  equation  for  U(x,t) 
takes  the  form  of  (1.51)  with  Un(t)  solved  in  (1.54).  Lastly, 
by  putting  x  «  1,  a  nonlinear  Vol terra  integral  equation  of 
the  second  kind  for  the  surface  temperature  U(l,t)  is 
obtained. 

As  in  the  previous  section,  the  integral  equation  for  the 
surface  temperature  will  be  explicitly  determined  for  two 
special  cases:  the  flat  plate  and  the  sphere.  Details  of  the 
derivation  of  the  solution  will  be  produced  in  the  case  of  the 
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flat  plate,  but  only  major  results  will  be  given  in  the  case 
of  the  sphere. 

gass  II  «i  =  1,  02  =  0,  03  »  -1,  h  =  1 

As  mentioned  in  section  1(C),  this  set  of  parameters 
corresponds  to  the  geometrical  configuration  of  a  flat  plate. 
Substituting  the  values  of  Sj,  ttj,  and  in  (1.49)  and  (1.50), 
Cl  equals  zero,  and  (1.50)  leads  to 

cosp^  =  PaSinp^  =>  -^  =  tanP^  ,  (1.55) 

where  cos  P„  ^  0. 

So,  the  family  of  orthogonal  eigenfunctions  are 

X^(x)  =  COS{PaX)  ,  (1.56) 

where  n  *=  1,2,3,...,  and  {P„)„.i“  is  the  set  of  distinct 
eigenvalues  which  are  the  roots  of  (1.55)  with  the  property 

0<Pi<P2<p3<- • • 


Next,  applying  the  finite  Fourier  integral  transform  of  the 
heat  equation  (1.1)  yields  (1.54)  in  terms  of  X„(x).  Using 
the  boundary  conditions 


dU(0,  t)  _ 

3x 


(1.57) 
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-  h  U*{l,t)  ,  (1.58) 

0  ,  (1.59) 

(1.60) 


and  the  fact  that 

Xl,'{x)  =  -  PX(x)  (1.61) 

produces  the  following  ordinary  differential  equation  for 
Un(t) 

+  PX(t)  ^  -  h  X^{l)U*(l,t)  .  (1.62) 

Note  that  (1.62)  is  a  first  order  linear  ordinary  differential 
equation.  We  find  the  solution  to  be 

U„{C)  =  -  h  dx  ,  (1.63) 

where 

U„{0)  «  (^g{x)Xg{x)  dx  .  (1.64) 

Jo 


X'(l)  +  X{1)  = 
X'(0)  =  0  , 
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Thus,  with  h  =  1,  the  integral  equation  for  U(x,t)  takes  the 
form 


=  e:.i 


^  -  (1)  ctv]x„{x) 


ix) 

Jo 


dx 


where  U„(0)  and  X„(x)  are  defined  by  (1.64)  and  (1.56), 
respectively.  Lastly,  by  putting  x  =  1,  the  integral  solution 
for  the  surface  temperature  U(l,t)  is  determined  to  be 


U(l,  t) 


( 1 )  p9(x)  X„  ix)  dx 

rx^(x)  dx 
Jo 


Jo - ^ - }  , 

f  xUx)  dx 
Jo 


(1.65) 


where  g(x)  is  the  initial  condition,  and  Xn(x),  and  are 
defined  as  above. 

Case  2;  =  0,  02  =  -1,  aj  =  1,  h  =  1 

In  this  case,  a  spherical  body  is  considered.  In  a 
similar  fashion,  the  family  of  orthogonal  eigenfunctions  can 
be  found  and  are  given  by 

X„[x)  =  sin(p^)  ,  (1.66) 
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where  n  »  1,2,3,...,  and  is  the  set  of  distinct  eigenvalues 
that  are  the  roots  of 

Pa  =  tan  (1.67) 


with  the  property 

0  <  Pi  <  Pj  <  Pj  <  ... 

After  applying  the  finite  Fourier  integral  transform  of  heat 
equation  (1.1)  with  respect  to  X„(x),  the  following  ordinary 
differential  equation  for  U„(t)  is  obtained 

*  PX(t:)  =  -  .  (1.68) 

Thus,  the  solution  of  equation  (1.68)  is 

£7a(t)  =  UjO)e'^‘^  -  h  dx  ,  (1.69) 


where 

=  f^g(x)Xg(x)dx  . 
Jo 


(1.70) 


So  therefore,  with  h  -  1,  the  integral  equation  for  U(x,t) 
takes  the  form 


,{x) 


[Oa(0)e'^“^-  (1)  Xjl) 

=  E;.a  ^ - TTT— - ’ 

Jfa  (X)  dx 
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where  U„(0)  and  X„(x)  are  defined  by  (1.70)  and  (1.66), 
respectively.  Lastly,  by  putting  x  =  1,  the  integral  equation 
for  the  surface  temperature  U(l,t)  becomes 


U{1,  t) 


(l)f^g{x)Xjx)  dx 

f^X^(x)  dx 
Jo 


Jo _ 

f^Xa(x)  dx 
Jo 


I 


where  g(x)  is  the  initial  condition,  and  X„(x)  and  are 
defined  as  above. 

E.  REMARKS 

The  solution  presented  above  is  not  complete  in  the  sense 
that  the  surface  temperature  is  only  determined  for  two  cases. 
The  solution  for  other  geometrical  configurations  can  be  found 
in  some  of  the  literature  listed  in  the  references, 
specifically  3,  5,  6,  and  11. 

The  surface  temperature  solutions  which  have  been  derived 
above  by  both  methods  fall  into  the  form 

U{l,t)  =4>(t)  -  +  £;;,^b*e-‘^*‘‘-’’}F[£/(l,T)]  eft,  (1.72) 

where  F  is  a  nonlinear  function  of  U(l,t),  and  c,,,  b^,  a,  and 
h  are  some  constants.  Equation  (1.72)  is  a  nonlinear 
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Volterra  integral  equation  of  the  second  kind.  4>(t)  is  a 
function  which  is  usually  called  the  "lag"  part  of  the 
integral  equation.  The  integral  in  (1.72)  is  often  referred 
to  as  the  "Volterra"  part  of  the  integral  equation.  In 
addition,  the  piece  within  the  braces  of  the  Volterra  part  is 
called  the  "kernel"  of  the  integral  equation.  As  these 
integral  equations  are  being  examined,  several  facts  about 
(1.72)  are  summarized  as  follows: 

1)  .  All  of  these  integral  equations  are  singular  because  as 
T  approaches  t,  the  kernel  blows  up  to  infinity. 

2)  .  All  of  the  infinite  series  satisfy  the  following 
property: 

If  f(t-T)  is  used  to  denote  an  infinite  series, 
then  lim^.  .  f(t)  =  constant,  thus  remaining  finite. 

3)  .  The  lag  part,  ^(t)  ,  and  the  kernel  of  the  integral 
equations  are  determined  by  the  geometry  of  the  body 
considered. 

The  above  "facts"  are  concluded  from  the  two  special  cases 
without  loss  of  generality.  In  each  of  the  next  three 
chapters,  a  different  numerical  method  for  solving  the  problem 
stated  in  section  A  will  be  introduced.  Both  the  method  of 
successive  approximations  and  the  Runge-Kutta  method  are 
numerical  techniques  used  to  deal  with  the  integral 
representation  of  the  problem,  whereas  the  finite  difference 
method  is  applied  directly  to  the  governing  equations. 
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IZ.  THE  METHOD  OF  SUCCESSIVE  APPROXIMATIONS 


A.  INTRODUCTION 

The  surface  temperature  of  a  body  subject  to  a  combined 
convective  and  radiative  boundary  condition,  as  seen  in 
Chapter  I,  is  given  by  the  solution  of  a  singular  nonlinear 
Vol terra  integral  equation  of  the  second  kind.  Since  the 
integral  equation  is  not  in  closed  form  and  is  nonlinear, 
numerical  techniques  seem  to  be  the  most  practical  way  to 
tackle  the  problem.  Over  the  past  twenty  years,  a  lot  of 
research  has  been  done  on  the  numerical  solution  of  an 
integral  equation  of  the  form 

:7(i,t)  *4>(t)  -  dx  .  (2.1) 

Among  the  existing  numerical  methods  for  solving  (2.1),  the 
method  of  successive  approximations  is  the  most  popular  one 
( see  [ 1 ] ) .  It  is  based  on  the  idea  that  the  set  of  successive 
functions  defined  by 

-  hJ^^k(t-T)F(ya(T) )  eft  ,  (2.2) 

where  k(t-T)  is  equal  to  the  term  in  braces  in  (2.1), 
converges  to  a  solution  of  (2.1)  in  every  finite  interval  of 
time.  In  the  following  section,  the  solution  method  will  be 
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outlined,  and  at  the  end  of  the  chapter,  general  comments  will 
be  made  on  the  technique. 

B.  OUTLINE  OF  THE  METHOD 

Consider  the  time  domain  in  which  integral  ec[uation  (2.1) 
is  to  be  solved.  Suppose  the  domain  is  partitioned  into  N 
inteirvals.  For  the  first  time  interval,  Ostiti  »  the 
approximate  solution  of  the  integral  equation  can  be  obtained 
by  using  the  iteration  procedure 

=^{t)  -  hf*'k{t-x)F{U„(l,x))  dx  (2.3) 

until  the  error  between  two  approximations  is  less  than  a 
predefined  small  number.  Next,  consider  the  second  time 
interval,  ti^titj.  In  this  interval, 

=^(e)  -  hf*'k{t-x)F(U„{l,x))  dx  (2.4) 

Jo 


can  be  broken  into 

=4>{t)  -  hf^^^k{C-x)F{U„{l,x))  dt 

-  hj^k[t-x)F{U„{l,x))  dx  .  (2.5) 

Since  U(l,t)  is  determined  for  O^tstj,  the  first  integral, 

hj^k{t-x)F{U„{l,x))  dx  , 
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is  known  (approximately) ,  and  thus  the  iteration  procedure  is 
only  needed  for  the  second  integral.  In  generrl,  for  the  i-th 
time  interval,  the  iteration  procedure  is  given  by 


=4>(t)  -  Jif‘''"k(C-x)F(U^(l,x))  dx 
Jo 

-  iif  ic(t-t)F(C7„(l,T) )  dt  ,  (2.6) 

•'  ^1-1 


where  t^.i^tsti. 

As  the  procedure  continues,  the  surface  temperature  is 
found  for  all  desired  times.  One  may  notice  that  as  the 
algorithm  is  carried  out,  the  singularity  of  the  Volterra  part 
of  the  integral  equation  creates  difficulty.  Appropriately, 
one  has  to  know  the  nature  of  the  singularity  which  the  kernel 
possesses.  To  illustrate  the  idea,  consider  the  integral 

dz  .  (2.7) 

J»  •Jh-z 

This  integral  is  often  found  in  the  integral  representation  of 
the  stated  problem.  Integral  (2.7)  possesses  a  sincfularity 
which  can  be  removed  by  the  use  of  the  transformation 

z  =  a  +  (b-a)  (l-X*)  .  (2.8) 
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Then,  by  using  a  suitable  Gaussian  quadrature  formula,  the 
integral  can  be  evaluated  accurately.  Normally,  one  usually 
comes  up  with  an  integral  with  a  stronger  singularity. 

The  iteration  procedure  outlined  above  needs  a  starring 
value.  Generally,  the  algorithm  will  converge  faster  to  the 
exact  solution  if  the  starting  value  is  close  to  the  exact 
solution.  Thus,  the  choice  of  initial  approximation  is 
crucial  for  convergence.  Based  on  the  fact  that  the  solution 
of  the  stated  problem  is  continuous,  one  can  choose  the 
temperature  at  the  previous  time  level  as  the  first 
approximation  of  the  method  when  a  small  time  step  is  used. 

The  method  of  successive  approximations  has  been  applied 
(in  Chapter  II)  to  solve  integral  equation  (2.1).  In 
particular,  a  method  used  to  tackle  a  simple  type  of 
singularity,  which  one  may  encounter  when  evaluating  the 
Volterra  part  numerically,  has  been  discussed.  Since  numerical 
integration  is  one  of  the  key  steps  in  the  method,  the  choice 
of  the  numerical  integration  scheme  does  affect  the  overall 
performance  of  the  algorithm.  One  can  improve  the  accuracy  of 
the  successive  approximations  method  by  appropriately  choosing 
a  numerical  quadrature  that  can  best  deal  with  the  singularity 
found  in  the  integral  equation.  Even  though  the  procedure 
outlined  above  may  seem  simple,  it  has  been  shown  that  the 
method  is  impractical  for  large  times  [ 3 ] . 
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III.  THE  RUNGS-KUTTA  METHOD 


A.  IMTRODUCTIOM 

This  chapter  considers  another  way  to  deal  with  the 
integral  equation  for  the  stated  heat  conduction  problem, 
namely  the  Runge-Kutta  method  which  was  first  introduced  by 
Crosbie  and  Viskanta  [5]  in  1968.  The  basic  idea  of  the 
method  is  based  on  an  approximation  of  the  kernel  by  a 
separable  kernel.  The  integral  equation  is  differentiated 
with  respect  to  time  and  transformed  into  a  nonlinear 
differential  equation.  The  Runge-Kutta  method  is  a  well  known 
numerical  scheme  for  solutions  of  ordinary  differential 
equations.  In  order  to  employ  the  method  the  surface 
temperature  at  a  desired  time  must  be  determined.  The  order 
of  approximation  of  the  method  is  determined  by  the  order  of 
the  ordinary  differential  equation.  What  differentiates  this 
method  from  the  other  numerical  schemes  is  instead  of  solving 
an  integral  equation  directly,  the  Volterra  integral  equation 
is  first  reduced  to  a  system  of  nonlinear  ordinary 
differential  equations  and  then  solved  numerically.  The 
method  is  not  exact  since  the  approximation  of  the  kernel  is 
not  practical  if  time  steps  are  small.  The  accuracy  of  the 
approximation  of  the  kernel  increases  with  time  and  order.  In 
the  next  section,  the  method  will  be  outlined  in  detail  as  it 
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is  applied  to  the  Integral  equation  (1.72) .  In  addition,  as 
an  example,  the  formulas  for  the  third  and  the  fifth  order 
versions  of  the  method  will  be  presented  explicitly. 

B.  OUTLINE  OF  THE  METHOD 

Consider  the  integral  equations  derived  in  Chapter  I . 
Generally,  the  integral  representation  for  the  dimensionless 
surface  temperature,  U(l,t),  of  the  body  that  we  have 
considered  can  be  written  as 

U{1,C)  =<|>(t)  -  f‘k{t-x)F{U{l,x))  dx  ,  (3.1) 

Jo 


where 

=  Po  .  i3.2) 

The  function  F(U(l,t))  is  the  surface  heat  flux;  the  o^'s  and 
Pk's  are  eigenvalues  and  coefficients,  respectively.  As  shown 
in  chapter  1,  the  infinite  series  k(t-T)  has  the  property 

lim^^JcC t)  =  Po  •  (3.3) 
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This  is  a  necessary  condition  for  an  integral  equation  to 
which  the  method  is  applied.  Now,  assume  ^(t)  is  a  bounded 
differentiable  function.  The  N— -order  approximation  of  k(t-T) 
is  given  by  taking  the  first  N  terms  of  the  infinite  sum.  So, 
(3.2)  becomes 

ic(t-t)  «  Po  +  5^^  .  (3.4) 


Substitute  (3.4)  in  (3.1)  and  let 

J*(t)  =  V*’F(£7(1,t)  )  dv  .  (3.5) 


U(l, t)  becomes 

t;(i,t)  =4>{t)  -  poj‘^V(y(i,t))  eft  - p*jr*(t)  .  (3.6) 

Differentiating  with  respect  to  time,  equation  (3.6)  becomes 

D<^Ml,t)  =4>‘^’(t)  -PoF(D(l,t)) 

*(F(D(1,  t) )  -  ajj*(t)]  .  (3.7) 

D<2)(l,t)  -  PoF‘iMc?{i,  t)  ) 

-  5^.1  P*[F‘i>  (£^(1,0)  -«*F(D(l,t))  +ajj*(t)]  ,  (3.7a) 
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*  {-i)-ai"j*{t)] .  (3.8) 

In  general,  the  N— -order  approximation  of  the  surface 
temperature  U(l,t)  is  determined  by  assuming  that  k(t-T)  in 
(3.1)  takes  the  form  of  (3.4),  in  which  only  the  first  N  terms 
of  the  infinite  sum  are  considered,  so  that  U(l,t)  is 
approximated  by  (3.6).  Then,  by  performing  N+1 

differentiations  of  (3.6),  the  resulting  system  of 
integrodifferential  equations  obtained  by  substitution  of 
1,  ...,  N+1  for  m  in  (3.8)  is  found  to  be 

-  PoF(cr(i,  0) 

^[F{U{1,  t))  -  ajj*(t)]  ,  (3.9) 

cr«Mi,t)  =<|>‘*Mt)  -  ^^F^^Hu{i,t)) 

-  J^iPjc[i^<^Mc/(i,  0)  -  alF(t?(i,  t))  +aij*(t)]  .  (3.10) 
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(1.  t)  -  (t)  -  t) ) 

-  rtx  picCi;  (-i)^ai^F<^-^-i)  mi.  t) ) 

+  (-l)*^aj^Jjt(t)l  ,  (3.11) 


£r<"*^>  (1,  t)  =  (t)  -  PoF^*”  iu{l,  t) ) 

-  Ei-i  (-D^oi^F'*-^'  (£7(1,  t) ) 

♦  (-l)»‘»«“*7lj(t)]  .  (3.12) 

To  eliminate  integrals  Ik(t),  where  ksl,...,N,  from  (3.12), 
consider  the  first  N  derivatives  of  the  surface  temperature 
which  are  given  by  (3.9)-(3.11) .  Rearrangement  of  the  terms 
in  (3.9)-(3.11)  yields 

J^^ajtp*J*{t)  ^U^^Ul,t)  -4><i>(t)  +  PoF(t7(l,  t)) 

+  5]^iP*F(rJ{l,  O)  ,  (3.13) 

-  J^.,atP*J*(t)  =C7<»J(l,t)  +  PoF<i)(c^(l,t)) 

+  5^iP*(F‘i»(Cf(l,t))  -  a*F(t^(l,  t))]  ,  (3.14) 


37 


+  (-l)V*^V<"-*-^»(Cf(l,t))]  ,  (3.15) 


(-1)"*"  5Xia*^pjrrjc(t)  =£7<"(i,t)  -♦‘'"(t)  +  Po^‘"‘"’(c^(i,t)) 
+  EJ.,  Pi[5]?li  (17(1.  t))]  .  (3.16) 


In  matrix  representation  we  obtain. 


Pi«i  P2*l 

-Pl«l  "P2«2  - 

*  •  • 

li(t) 

liit) 

s 

p  « 

t  s  • 

(-l)^P,ar*  (-i)"P2«r’  ...  (-D^Pirar* 
(-i)"*iPiaf  (-D^^'Pjar  ...  (-1)"*"PX* 

lif-i,  ( t) 

J*(t) 

^#-1 

where  Bi....,Bn  are  defined  as 

-  C7‘i>  (1,  t)  +  poF(t7(l,  t)) 

P*F(C7(l.t))  ,  (3.18) 
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*  J^^Pj,[Ff^Ucr(i,e))  -  a}F(cr(i,t))]  ,  (3.i9) 


Sh-i  =  <1/  t)  '  (t)  +  (t^d,  t) ) 

-  n.  P4S';:  <  -i)^a|^F<^-2-^Mcr(i,  t))]  ,  (3.20) 


*  n.1  pitKij  (-1) It'd- 1) )]  •  (3.21) 

Now,  let 

Pltt^  Pjtt^  ...  Pj/X^ 

“Pl*l  ”P2*2  '•• 

:  :'  ...  5  .  (3.22) 

(-i)"P,ar"  (-i)^P2«r"  ...  (-l)'^P^ar* 

(-D^^ip^a?''  (-1)^*%«I'' ...  (-l)''*"Pi^f 

By  Cramer's  rule,  Ii(t) , . . .  ,lN(t)  can  be  expressed  as  a 
quotient  of  N  x  N  determinants,  given  by 
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(3.26 


where 


Det  {A) 


Pia? 

p2«i 

-Piat 

-P2«2 

”Pi»®W 

(-D^Pia^" 

•  ••• 

{-i)"p^r" 

(-l)^*^Piaf 

(3.27) 


Next,  by  using  only  the  fiindamental  properties  of 
determinants,  (3. 23) -(3. 26)  can  be  simplified  as  a  quotient  of 
(N-l)x(N-l)  determinants  (see  appendix-A  for  further  detail)  . 
Define 


Det  {A')  = 


-(«^a?) 

(aS-at) 


-(«3-«?) 

(-i)*f(af^-^-er’*) 

(-l)'*Mar*-ar*) 


-(al-al) 

(«*-«!) 


(3.28) 


41 


The  formulas  are  found  to  be 

Ji(t)  - 

(ajr-aa) 

PittfDetCAO 

(3.29) 

and  for  k  *=  2,  ...,  N-1, 

(continued  on  next  page) 


(Bj+ttzBi)  -(ai-al) 

«  • 

•  •  ••• 

{B^.i-(-i)"ai^'^Bi)  (-i)^(gf^^-g|"‘^  ... 

(B„-{-l)"*"a|"-"BJ  (-l)^*i{ar‘*-gf-")  ... 
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(3.30) 


-(ai-ow-i)  ...  -ial-i-ol-i) 

(tti-ajlr-i)  ...  (•Bj'Oft-i'Si) 

(-l)^*Mtti"~^-ay-I^)  ...  {Bt,-{-i)^^ci^l^B^) 

^^sDetiA') 

(3.31) 

Thus,  Ik(t),  where  k  -  1,  ...,  N,  in  equation  (3.12)  can  be 
determined  explicitly  by  formulas  (3.29)  through  (3.31). 
Hence,  the  integrodifferential  equation  (3.12)  is  reduced  to 
a  (N+1)—  order  nonlinear  ordinary  differential  equation 

-  po^‘«  (cr(i,  O) 

■  TLi  ( -!)"«"'''*■'’  («i.  t) ) 

♦  ,  (3.32) 

with  the  initial  values 

a(l,0)  =  ♦(O)  ,  (3.33) 

C?<‘Ml,0)  =4><^M0)  -  ^^FiUil.O)) 

^t[F{U{l,0))]  ,  (3.34) 
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-  2!.i  (tr{l,0))  -  O*F(Cr(l,0))]  ,  (3.35) 


-  5^^  (C?(1,0))]  ,  (3.36) 

Which  can  be  obtained  by  putting  t  »  0  in  (3.9)  through 
(3.12),  andli(t),  ...,  InCt)  are  determined  by  formulas  (3.29) 
through  (3.31). 

The  Runge-Kutta  method  is  then  applied  to  the  nonlinear 
ordinary  differential  equation  (3.32)  with  initial  conditions 
(3.33)  through  (3.36).  Hence,  the  Nil;  order  approximation  of 
the  surface  temperature  will  be  given  by  the  numerical 
solution  of  a  (N+1)—  order  nonlinear  ordinary  differential 
equation. 

As  an  example,  in  the  next  section,  formulas  for  the  third 
and  the  fifth  order  approximation  of  the  method  will  be 
presented.  Details  of  the  derivation  of  the  equations  will 
not  be  produced,  and  only  major  results  will  be  given. 
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C.  THE  THIRD  ORDER  APPROXIMRTXOM 


The  third  order  approximation  of  the  surface  temperature 
of  (3.1)  is  given  by  the  numerical  solution  of  the  following 
fourth  order  nonlinear  ordinary  differential  equation 

=  <!><*>(  t)  -  (Po  +  Pi  +  Pa  +  P3)F<3J(C^(l,t)) 

+  (Pia?  +  Pa®!  +  (:^(1,  t)) 

-  (Piflfl  +  *  p3a5)F‘iMt^(i,  t)) 

+  (Piflf?  +  p2«l  +  P3a|)F(u(i,  t)) 

-  {Pia?-ri(t)  +  Paaf-TaCt)  +P3a|J3(t))  (3.37) 

with  initial  conditions 

£7(1,0)  =  ^(0)  ,  (3.38a) 

£7<i’(l,0)  =<|>‘i>(0)  -  (Po  +  Pi  +  Pa  +  P3)i^(<l>(0))  ,  (3.381?) 

C7«2)(i,o)  =  <|>‘2)(o)  -  (po  +  Pi  +  Pa  +  P3)F<i)«|>(0)) 

+  (Pitti  +  P2«l  +  P3«|)i^(<>(0)  )  ,  (3.39) 


and 
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J7<3)(1,0)  -  (Po  +  Pi  ♦  P2  + 


Il(t), 


where 

B, 


+  (PiOi  +  P2«l  +  P3a3)F‘^’ (♦(0) ) 

-  (Plttt  +  P2«2  +  P3«3)i^(4>{0))  .  (3.40) 


EjCt),  and  l3(t)  in  (3.37)  are 


Ji(t) 


(Bj+alBi)  -(a3-«2) 
(Bj-aaBi)  (a3-«2) 


Pitti  Det{A') 


(3.41) 


I,(t) 


(a^-af)  (Ba+a^Bi) 
(at-a5)  (Bj-a^Bi) 


Patti  Bet  (A') 


(3.42) 


lAt)  =  - 


(tt|-tt|)  (Bj+ttlBj) 

(ttl-tt|)  (S3-tt|Bi) 


PjttI  Dec  {A') 


(3.43) 


=  B<i>(l,t)  -  <|)‘‘MC)  +  (Po  +  Pi  +  P2  +  p3)F(B(l,t))  , 
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-  (Po  +  Pi  +  Pa  +  P3)F<i>(C?(l,t)) 


+  (Pitt?  +  P2«2  ♦  P3a|)F{:7(i,  t))  , 


B3  =  cr»)(i,t)  +  (Po  +  Pi  +  Pa  +  P3)F<"’(tr(l,t)) 

-  (Pi«i  +  p2ai  +  hal)F^^Hu{i,t)) 

*  (Piai  +  P2a2  +  p3aS)F(cr(i,  t))  , 

and 

-(a|-ai)  -(af-a^ 

DetUO  =  .  4  • 

iat-ctx)  (tta-ai) 


O.  THE  FIFTH  ORDER  RFFROZZMATXON 

The  fifth  order  approximation  of  the  surface  temperature 
of  (3.1)  is  given  by  the  numerical  solution  of  the  following 
sixth  order  nonlinear  differential  equation 
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-  <!»<«’(  t)  -  (Po  +  Pi  +  p2  ♦  p3  +  P4  +  P5)i^‘®MCr{l,t)) 

+  (Pi*!  +  P2«2  +  p3«i  +  P4a4  +  p5ai)i^‘*’  mi,  t)  ) 

-  (PiOfl  +  P2«2  *  P3«3  +  p4«4  +  p5a5)F‘3)  {U{1,  t)  ) 

+  (Piat  +  P2a2  +  Patti  +  P4tt|  +  p5tt|)F‘2)  {u{i,t)) 

-  (Pittl  +  Patt®  +  Pattf  +  P4a|  +  P5af)F<i»  (C^{1/ 1) ) 

+  (Piai°  +  pjttl®  +  Patti®  +  P4tt|°  +  p5tt|°)F(cr{l,  t) ) 

-  (Pitti"ji(t)  +p2tt|^J2{t)  +P3tt|^J3(t) 

+  P4tt“j4(t)  +P5tt|^J5(t))  ,  (3.44) 

with  Initial  conditions 

mi,Q)  =  4>(0)  , 

C7<i>(l,0)  =<|><i>(0)  -  (Po  +  Pi  +  P2  +  P3  +  P4  +  P5)F(<1>(0))  , 

Cr‘2)(i,0)  =<|><2){0)  -  (Po  +  Px  +  Pa  +  P3  +  P4  +  P5)^’<'M4>(0)) 
+  (Pitt?  +  Pjttl  +  Patti  +  p4tt|  +  Pgttl)  F(4>{0)  )  , 
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C^<3>(1,0)  =  <|><3)(t)  -  (Po  +  Pi  +  P2  +  P3  +  P4  +  P5)F‘*)(<K0)) 

+  (Pittf  +  P^af  +  Pattf  +  P^aJ  +  P5a|)F<iJ  (♦(0) ) 

-  (Pitti  +  pjtta  +  Patt^  +  P4al  +  p5a5)^’(4>{0) )  / 

U^*^{l,Q)  =<!»«<>  (0)  -  (Po  +  Pi  +  p2  +  Pa  +  P4  +  P5)^‘'M<|)(0)) 

+  (pitti  +  p2ai  +  Pattf  +  P4aJ  +  p5a|)F<2>  (<K0) ) 

-  (Pittl  +  P2a|  +  Patta  +  P4«J  +  Pgat)  F<^>  (4>(0) ) 

+  (Pittt  +  p2a|  +  Pattt  +  P40$  +  p5at)F(4>(0) )  , 

c;<5)(i,0)  =  4><®‘(0)  -  (Po  +  Pi  +  p2  +  Pa  ■*■  p4  +  p5)i=’^^’(4>(0)) 

+  (Pitti  +  pjal  +  Pattf  +  p4a4  +  Psal)^^^)  {4>(o)) 

-  (Pittt  +  P2«2  +  Patt^  +  p4aj  +  P5a5)F‘2’ (<1>(0) ) 

+  (Piat  +  Patti  +  Patti  +  P4ttt  +  P5ttt)F‘^’ (4>(0) ) 

-  (Pitt?  +  Patt?  +  Patti  +  P4a?  +  P5o|)F(4>(0)  )  . 

i(t),  l2(t),  l3(t),  I<(t),  and  l5(t)  in  3.44  are 
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Jl(t)  * 


J3(t) 


It(C) 


(Bj+a^Bi) 

-(a3-o|) 

“(aj-«i) 

-ial-al) 

iat-at) 

(ol-aj) 

iat-at) 

(B^+afBi) 

-(a|-aS) 

-iat-al) 

-iat-at) 

(Bs-alBj) 

ial-al) 

ial-al) 

(ai-o|) 

Pitti  Det{A') 

-(a?-al) 

(Bj+a^Bj 

“(a^a^) 

-(oti-os) 

(al-a^) 

(Bj-a^Bj) 

(al-a5) 

iat-al) 

-(«!-a5) 

(B^+a^Bi) 

-(aS-a^) 

-ial-at) 

(Bj-a^Bi) 

(aj-af) 

Pjtta  Det{A') 

-(ai-a*) 

-(«|-«4) 

(Bj+a^Bj) 

-(«|-a^) 

iai-at) 

(a^«^ 

(«^«J) 

-(al-at) 

-(o|-at) 

(B^+aSBj) 

-(o|-aS) 

(aJ-aS) 

(ol-aj) 

(Bj-aSBi) 

(a|-a5) 

Pjtta  DetiA') 

-(a!-ai) 

-ial-al) 

-(a|-ai) 

(Bj+aiBj) 

(al-a?) 

{al-at) 

{al-at) 

(B3-05B1) 

-(a^-af) 

-(al-a|) 

-(af-af) 

(B4+0IB1) 

(a$-a|) 

(al-al) 

(a;-a|) 

(Bj-agBi) 

net  {A') 


(3.45) 


(3.46) 


(3.47) 


(3.48) 
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I^it) 


-(ai-a*)  -(af-aj)  -(of-oj)  (^+04^1) 
(aj-oj)  (tta-aj)  (aj-aj)  (Bj-oJBi) 
-(af-at)  -(a2-«$)  -(a|-a$)  {B^+alB^) 
(ttS-a!)  (a|-ttS)  (aS-tt")  {B^-alB^) 
Psttl  Det(A') 


(3.49) 


where 


Dot  {a') 


-(«!-«?)  -(al-al)  -(al-al)  -{al-al) 
(tta-al)  (tta-al)  (aj-at)  {at-a{) 
-(a|-al)  -(«!'«$)  -{a$-at)  -(a|-af) 
(a®-a?)  (af-«?)  (oj-a?)  (a|-a?) 


(1,  t)  -  (t)  +  (Po  +  Pi  +  P2  +  P3  +  P4  +  p5)i^(£^(i/ 1) ) 

Bj  =  +  (Po  +  Pi  +  P2  +  p3  +  P4  +  P5)i^‘"’ (a(i,  t) 

-  (Pl«l  +  P2«l  +  P3«l  +  P4«4  +  Ps^l)F{U{l,  t))  , 

B3  =  C7<3)  (1,  t)  -  (t)  +  (Po  +  Pi  +  P2  +  P3  +  P4  +  Ps)^**’  t) 

-  (pl«l  +  P2«2  *  P3«3  +  P4«4  +  PjOi)  t)  ) 

+  (Pittt  +  Patta  +  Patt^  +  p4«4  +  P5O5)  F(£r(l,  t) )  , 
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(1,  t)  -  +  (Po  +  Pi  +  Pa  ♦  P3  ♦  p4  +  p5)^‘^’  (^<1/  t) 

-  (Pitti  +  Patt*  +  Patti  +  p4tti  +  p5tt|)F‘*>  (1/(1/ t)  ) 

+  (Pjtti  +  Patti  +  Patti  +  P^ttJ  +  P5tti)F‘^>  (£/(!/ t)  ) 

-  (Pitt!  +  Patti  +  Patti  +  p4tt|  +  p5tt|)F(C7(l,  t))  , 

and 

Bg  =  (1,  t)  -  (t)  +  (Po  +  Pi  +  Pa  +  Pa  +  P4  +  (U{1,  t) 

-  (Pitti  +  Patti  +  Patti  +  P4tti  +  P5tti)F‘3M:7(l,  t)  ) 

+  (Pitt!  +  Patti  +  Patti  +  P4ttJ  +  PgttI)  t)  ) 

-  (PiOf!  P2«l  +  P3«t  +  P4a!  +  P5tt|)F‘^)(t^(l,  t)) 

+  (Pi«!  +  Pa*!  ♦  Pa®!  +  p4aS  p5tti)F(C7(l,  t))  . 


E.  REMARKS 

In  this  work  the  Runge-Kutta  method  is  applied  to  solve 
the  integral  equation  resulting  from  the  heat  conduction 
problem  with  combined  convection  and  radiation.  In 
particular,  the  nonlinear  ordinary  differential  equation  has 
been  determined  for  both  the  third  and  the  fifth  order 
approximations.  It  may  be  observed  that  this  method  is  not 
very  practical  for  calculating  the  temperature  at  small  time 
steps.  The  reason  is  that  the  smaller  the  time  one  takes,  the 
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more  terms  will  be  represented  correctly,  which  in  turn  may 
result  in  a  high-order  nonlinear  ordinary  differential 
equation  with  a  very  large  number  of  terms.  The  number  of 
terms  could  grow  to  infinity.  Thus,  the  method  is  usually  used 
to  compute  the  surface  temperature  at  large  times  where  the 
temperature  distribution  is  in  a  steady  state. 

In  the  following  chapter,  we  will  take  another  approach 
using  a  different  numerical  method,  namely  the  finite 
difference  method.  This  method  is  different  from  the  previous 
numerical  techniques  in  that  instead  of  solving  the  integral 
equation,  it  approximates  the  partial  differential  equation 
and  the  boundary  conditions  directly. 
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IV.  THE  TIHITB  DimSEIICB  METHOD 


A.  ZMTBODUCTZOH 

The  basic  idea  of  the  finite  difference  method  is  to 
transform  a  continuous  model  into  a  discrete  system  by 
replacing  the  continuous  domain  in  the  model  with  a 
den\imered}le  domain.  In  applying  this  idea  to  differential 
equations,  all  the  derivatives  in  the  equation  are  simply 
replaced  by  finite  difference  approximations.  Thus,  the 
unknowns  in  the  difference  equation  have  a  countable  domain, 
and  the  resulting  discrete  system  is  solved  numerically. 

In  the  theory  of  numerical  analysis,  the  significance  of 
the  computed  solution  of  a  finite  difference  scheme  in 
relation  to  approximating  the  exact  solution  depends  upon 
three  elements.  They  are  consistency,  convergence,  and 
stability.  Consistency  is  a  condition  used  to  assure  that  as 
AX  (the  spacing)  approaches  zero,  the  truncation  error  of  the 
scheme  also  goes  to  zero.  It  implies  that  the  finite 
difference  can  be  an  arbitrarily  accurate  approximation  to  the 
derivative.  Convergence  of  the  approximation  assures  that  if 
AX  goes  to  zero,  the  difference  between  the  computed  and  the 
exact  values  also  goes  to  zero.  In  other  words,  any  desired 
accuracy  of  the  approximated  solution  can  be  achieved.  The 
last  element  is  the  stability.  The  stability  of  a  scheme 
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concerns  the  growth  of  the  errors  found  in  the  calculations 
which  are  needed  to  solve  the  system  of  linear  equations.  A 
scheme  is  said  to  be  conditionally  stable  if  the  roundoff 
error  does  not  amplify  if  the  time  step  is  under  a  critical 
value  which  is  determined  by  the  differential  equation 
considered.  In  the  Lax  Equivalence  Theorem,  the  relationships 
of  these  three  conditions  are  stated.  It  says  that  given  a 
properly  posed  initial  value  problem  and  a  finite  difference 
scheme  which  satisfies  the  consistency  conditions,  stability 
is  the  necessary  and  sufficient  condition  for  convergence. 

There  are  many  difference  approximations  and  methods  for 
solving  discrete  systems  that  are  available  in  numerical 
analysis.  Different  choices  of  approximation  and  methods  of 
solving  the  system  will  lead  to  differing  degrees  of  accuracy 
in  the  approximation  of  the  solution.  This  chapter  will  only 
focus  on  a  particular  finite  difference  scheme  used  to 
approximate  the  governing  partial  differential  equation  in  the 
stated  problem  and  an  algorithm  for  solving  the  discretized 
system. 

B.  CRAMX-MICHOLSOM  SCHEME 

Suppose  a  lies  between  Xq  and  X{  and  t  z  to,  where  Xq  and 
Xf  are  some  initial  and  final  x-coordinate  which  brackets  the 
location  of  concern.  Let  ax  and  At  be  increments  of  x  and  t, 
respectively.  The  x-t  space  can  be  partitioned  into  a  grid 
network  in  which  the  points  are  given  by  x  =  Xq  +  jAX  and  t  = 
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to  +  nAt,  where  j  >  0,1,2,...,N,  with  N  being  the  nuiaber  of 
nodes,  and  n  «  0,1,2  ...  .  When  ax  and  At  are  constants,  the 
mesh  obtained  is  uniform,  and  the  temperature  at  x  »  Xq  + 
iAX,  written  as  X;,,  and  t  «  to  +  nAX,  written  as  t„,  is  denoted 
by  U^". 

As  previously  mentioned,  there  are  several  ways  of 
choosing  a  finite  difference  operator  for  replacing  the 
derivatives.  If  the  average  of  the  forward  and  backward 
difference  schemes  is  used  for  the  space  discretization  and 
the  forward  difference  scheme  is  written  about  the  point  Xi, 
t^.H,  the  governing  partial  differential  equation  becomes  a 
second  order  accurate  (in  both  x  and  t)  finite  difference 
equation.  It  is  given  by 

+  (-2  -  2p)£^"^  +  =  -  Uf-i  +  (2  -  2P)17/  -  Ujit  ,  (4.1) 

where 


(which  is  the  well  known  Crank-Nicholson  scheme) . 

Since  it  is  of  second  order,  the  truncation  error 
associated  with  (4.1)  is  on  the  order  of  o(ax^  +  At^)  .  Notice 
that  the  temperature  at  time  t„.i  is  a  function  of  unknown  and 
known  temperatures  at  six  of  the  ten  points  shown  on  the  Fig. 
4.1. 
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(4.2) 


N'l  foz  case  l 
N-i  foz  case  2 


To  ensure  that  the  oscillation  is  eliminated,  the  implicit 
backward  finite  difference  scheme  (which  is  satisfactory  with 
all  types  of  boundary  conditions)  is  adopted  at  the  boundary, 
X  >  1.  The  equation  at  x  «  1  is  given  by 

US:}  +  (-2  -  *  us:^;  =  -  ,  (4.3) 

where  P  is  as  before. 

There  is  a  fictitious  point  outside  the  computational  domain 
in  (4.3),  that  is,  the  unknown  temperature  at  N+1  is  denoted 
as  To  eliminate  that  point,  use  a  difference  method  to 

approximate  the  derivative  in  the  radiative  boundary  condition 
(1.2c)  because 

TfO*l  _  TfO*! 

-"2Tx—  ■ 


where  F  is  the  right  hand  side  of  (1.2c).  Algebraically 
manipulating  (4.4)  yields  the  following  equation 

uS*i  =  US-i  *  2Axa3US*^  +  2axF{uS*^)  .  (4.5) 

Substituting  (4.5)  into  (4.3),  the  resulting  expression 
becomes 
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20|!?-1  +  (-2  -  p  +  2AXa^)uS*^  =  -  pCTjf  -  2axF{US*^) 


(4.6) 


which  is  a  nonlinear  equation  in 

Observe  that  (4.1)  and  (4.6)  constitute  a  set  of 
simultaneous  equations  at  each  time  step.  In  matrix 
representation,  the  resulting  system  is  of  the  form 

AU  ss  B  ,  (4.7) 

where  A  is  a  tridiagonal  matrix,  B  is  a  vector  of  all  the 
known  values  found  in  each  equation,  and  U  is  a  vector  of  the 
unknown  temperatures  at  each  space  node  at  a  particular  moment 
of  time.  So,  for  each  time  level,  the  transient  temperature 
is  given  by  the  solution  of  a  system  of  equations. 

The  Thomas  algorithm  can  be  used  to  solve  a  tridiagonal 
system  of  linear  equations.  Clearly,  all  the  equations  in 
(4.7)  are  linear  (except  the  last  one).  The  first  half  of  the 
algorithm,  as  given  in  appendix<-B,  can  be  directly  applied  to 
the  system  except  for  the  case  where  i  »  N.  In  the  case  of  i 
»  N,  substituting  d(N)  in  the  first  Do-loop  in  the  expression 
right  after  the  first  loop  yields 

,wj*i  _  d{lt)  -  ratio  *  d{N-l)  tA  q\ 

^ - F(W) -  ’  ' 
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which  implies 


-pCTj?  -  24xF(C^^)  -  ratio  ♦  cl(N-l) 
b{N) 


(4.9) 


with  b(N),  d(N~l),  and  ratio  computed  in  the  first  Do- 

loop  (reference  to  appendix-B) .  Now,  (4.9)  can  be  rewritten  as 


-P£7j?  -  2axF(uS*^)  -  ratio  ♦  d{N-l) 

Fun 


US*^  =  0  .  (4.10) 


Let  the  left  hand  side  of  (4.10)  be  represented  by  f.  It 
follows  that 

f{US*^)^0  .  (4.11) 

Thus,  the  update  of  the  surface  temperature  is  the  solution  of 
the  nonlinear  equation  (4.11). 

In  the  following  section,  the  cases  for  a  flat  plate  and 
a  sphere  will  be  considered  to  obtain  the  respective 
tridiagonal  systems. 

C.  TWO  SPECIAL  CASES 
1.  The  Flat  Plate 

The  parameters  corresponding  to  this  case  can  be  found 
in  Chapter  I.  Applying  the  finite  difference  method  outlined 
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above  to  the  governing  equations  leads  to  the  following 
results . 

Consider  the  Crank-Nicholson  scheme  for  j  »  0,  N-1. 

For  j  “  0, 

+  (-2  -  =  -  Ufi  +  (2-2P)C;f  -  .  (4.12) 


To  eliminate  the  fictitious  points,  the  boundary  condition  at 
X  ^  0  in  discretized  form  is  taken  to  be 


Ua*l  _  jjn*l 

2ax 


=  0 


(4.13) 


Thus, 

,  (4.14) 

and 

=  Ufi  .  (4.15) 

Substituting  (4.14)  and  (4.15)  in  (4.12)  produces 

(-2  -  2P)0'o*^  +  2C^*^  =  (2-2P)0f  -  2U°  .  (4.16) 
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For  j  »  1,  ,  N-1 

+  (-2  -  2P)C^*^  +  =  -!;“+  (2-2P)£^  -  t/j"  , 

US:2  +  (-2  -  2p)£7j?.Y  +  US*^  =  -  US.2  +  {2-2p)C;^.i  -  CTjJ'  . 

When  j  «  N,  as  shown  before,  the  equation  becomes 
2£)5?.*i  +  (-2  -  p  +  2lX(i^)U^*^  =  -  pCTjf  -  2axF{US*^)  . 


In  matrix  representation,  with  initial  values 


crj  »  1 

,  where 

j 

0 ,  .  .  .  ,  N  f 

we  have 

(-2-2p) 

2 

0 

••• 

0 

1 

(-2-2p) 

1 

0 

0 

0 

1 

(-2-2P) 

1 

0 

0 

0 

... 

... 

0 

A  = 

• 

! 

... 

... 

0 

; 

: 

... 

... 

...  ; 

0 

0 

... 

1 

(-2-2P)  1 

0 

0 

... 

0 

2  (-2-p+2AXa3) 

(4.17) 


(4.18) 


(4.19) 
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(2-2p)Uo"  -  2i7i” 

at' 

-Uo  +  (2-2P)U"  -  U2 

at' 

* 

,  and  U = 

• 

-U^.2  +  (2-2P)U;f.i  -  Uj? 

us:i 

-puj?  -  2j>jcF{uS*^) 

2 .  Thtt  Bphmrm 

The  parameters  again  can  be  obtained  in  Chapter  I  and 
will  not  be  repeated  here.  Using  the  finite  difference  method 
outlined  above  with  the  governing  equations  leads  to  the 
following: 

Consider  the  Crank-Nicholson  scheme  for  j  =  1,  N-1. 

For  j  =  1, 

+  (-2  -  ^  -  Uq  *  (2-2p)£^  -  U2  .  (4.21) 

However, 

=  0  and  Uq  =  0  .  (4.22) 

Substituting  (4.22)  in  (4.21)  produces 

(-2  -  2^)ur^  +  =  (2-2P)l^  -  U2  .  (4.23) 


For  j  =  2,  ...»  N-1 
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£^♦1  +  (-2  -  2^)US*^  +  +  {2-2P)l}^  -  Uf  , 


(4.24) 


US:2  +  (-2  -  2P)t7^:^  +  +  (2-2P)C^?.i  -  .  (4.25) 

When  j  «  N,  as  shown  before,  the  baclcvard  scheme  is  used. 
Thus, 

2USlt  +  (-2  -  p  +  2Axa^)  US*^  =  -  pUjf  -  2axF[uS*^)  . 

In  matrix  representation,  with  initial  values 
Uj  -  j/bjc  ,  where  j  =  1,  ...  ,  N  , 
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B 


and  U 


(4.37) 


{2-2^) 

-U^  +  {2-2f^)Ui  -  Uf 


UT^ 


-U^.2  *  {2-2P)£)r/.i  -  CTjf 
-  2axF{uS*^) 


D.  STABILITY 

Even  though  the  backward  analog  is  implemented  on  the 
radiative  boundary,  according  to  the  nvimerical  experiments, 
the  method  still  suffers  from  the  problem  of  oscillations  when 
a  large  time  step  is  imposed.  As  far  as  the  author  is  aware, 
not  a  single  formula  has  been  developed  for  the  stedsility 
criteria  of  an  implicit  finite  difference  scheme  with 
nonlinear  boundary  conditions.  However,  two  stability 
formulas  of  the  related  problems,  which  are  found  in  the 
literature  [7],  can  serve  as  a  guideline  in  choosing  the  time 
step  for  the  problems  considered.  The  first  one  is  due  to 
Lawson  and  Morris  [7].  They  deduce  the  stability  criterion 
for  the  Crank-Nicholson  equation  with  linear  boundary 
conditions  as 

At  <  .  (4.28) 

X 

Another  stability  criterion  is  due  to  Milton  and  Goss  [9]  who 
applied  the  laws  of  thermodynamics  in  developing  the  stability 
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requirement  for  an  explicit  finite  difference  scheme  with 
nonlinear  boundary  conditions.  It  turned  out  that  the  time 
step  required  for  the  stability  is  restricted  by 


At  2; 


(ax)  *inax{t7|f} 


(4.29) 


where  the  maximum  is  taken  over  all  n  and  where  Un°  can  be 
found  by  setting  the  following  function 

=  A  - 


equal  to  zero,  where 


4^,  =  ur  -  uS  , 


_  2US.^ 
(ax)* 


o  _  2Iax  -  1] 


2ax 

(ax)* 
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Because  this  formula  is  not  very  practical  in  actual  use, 
(4.28)  will  be  chosen  as  a  guideline  for  selecting  the  time 
step  of  the  method. 

E.  REFINEMENT  OF  PARTITION  AND  EXTRAPOLATION  TECHNIQUES 

The  partition  of  the  domain  covered  has  a  great  influence 
on  the  accuracy  of  the  solution  obtained.  The  choice  of  grid 
points  is  determined  by  knowledge  of  the  problem  and  by 
numerical  experimentation.  Here,  two  ways  of  improving  the 
accuracy  of  the  finite  difference  method  are  presented,  and 
one  of  the  two  is  chosen  to  be  implemented  in  the  numerical 
methods . 

One  way  to  improve  the  accuracy  is  the  so  called 
prolongation.  We  first  solve  the  problem  using  one  spacing 
and  then  refine  the  paxrtition  and  then  repeat  the  computation. 
If  the  comparison  shows  large  differences,  the  process  is 
repeated  for  smaller  and  smaller  grid  sizes  until  a  desired 
accuracy  is  achieved.  This  method  may  result  in  a  prolonged 
computational  time  for  the  solution. 

The  second  way  is  called  extrapolation.  The  simple 
ingenious  idea  of  the  technique,  which  dates  back  to 
Richardson  in  1910,  is  the  following: 

One  solves  the  same  type  of  problem  over  a  prescribed 
interval,  for  example  [0,1],  several  times  with  successively 
smaller  step  sizes.  Thus,  one  obtains  a  sequence  of 
approximations 
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Y  /  y  (i/iii)  /  . . . 

for  a  given  sequence  of  step  sizes 

Aq  >  Jjj  >  ...  >0  . 

The  successive  step  size  hj  is  often  defined  in  terms  of  an 
input  step  size  h  by 

hj  =  ,  i  *  0,1,2.  ...  .  (4.31) 

Thus,  any  step  size  sequence  (h^)  can  be  characterized  by  the 
associated  integer  sequence  {n^}.  The  following  are  some 
examples  of  Integer  sequences: 

{1,2,4,8,16,32,...}  (Romberg  sequence) 

{1,2,4,6,8,12,...}  ( Bui irsch  sequence) 

{1,2,3,4,...}  (harmonic  sequence) 

So,  the  numerical  solution  at  x  is  computed  for  a  sequence  of 
step  size  hi  and  denoted  by  Ti,o  *  Y(irl^i)>  Then,  the 
extrapolation  tableau. 

Too 

T20»T2i  ,T22 

is  calculated  for  x  according  to  two  types  of  commonly  used 
extrapolation  schemes 

a) .  Aitken-Neville  algorithm 
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For  1  B  1,2,...  and  k  *  l,2,3,...i 


(4.32) 


b) .  Rotational  extrapolation 

For  i  *  1,2,...  and  k  »  1,2,3,.. .,i 

li.-i  =  0  , 


Tj.it-i  -  r, 


i-i.jt-i 


-1 


(4,33) 


In  this  study,  extrapolation  scheme  (4.32)  with  the  Romberg 
sequence  will  be  implemented  when  the  finite  difference  method 
is  used  to  find  the  numerical  solution  of  the  stated  problems. 
It  should  be  noted  that  if  this  extrapolation  scheme  is  used 
the  computational  time  will  increase  exorbitantly. 


70 


V.  HDMBRXCAL  RESULTS 


A.  INTSOOUCTIOH 

The  problem  described  in  section  1(B)  was  solved 
numerically  for  two  special  cases,  namely,  the  flat  plate 
(ai=l,  a2«0,  a3=-l,  h=l)  and  the  sphere  (ai=0,  a2=-l,  «3=l, 
h«l) .  Since  a  lot  of  numerical  results  of  the  problem 
computed  by  successive  approximations  method  are  avail2d>le  in 
some  of  the  papers [3, 4],  in  this  thesis,  only  the  Runge-Kutta 
method  and  the  finite  difference  method  are  employed  to  the 
problem  for  study.  Programs  are  written  in  Fortran  77  using 
the  Amdahl  5990  model  500  mainframe  computer  and  are  set  up  to 
allow  input  for  the  time  step.  Thus  one  can  approximate  the 
maximum  time  step  that  can  be  used  in  a  particular  numerical 
method.  All  calculations  are  done  using  double  precision 
arithmetic  yielding  12 -digit  accuracy.  Numerical  results 
generated  by  the  methods  are  compared  and  discussed. 

The  Runge-Kutta  and  the  finite  difference  methods  are 
implemented  to  solve  both  special  cases.  In  particular,  three 
different  order  approximations  of  the  Runge-Kutta  method  are 
programmed  to  solve  the  integral  equations  derived  by  the 
Laplace  transform  method.  Inefficiency  of  a  high  order  Runge- 
Kutta  method  motivates  the  use  of  the  finite  difference 
technique.  Again,  the  method  is  implemented  in  both  cases  for 
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various  time  steps.  Some  of  the  nimerical  results  are 
tabulated  and  plotted  in  such  a  way  that  a  comparison  can  be 
made.  Notice  that  the  Runge-Kutta  method  is  not  applied  to 
the  integral  equations  obtained  by  the  eigenvalue  expansion 
method.  The  reason  is  the  lag  parts  of  those  integral 
equations  diverge  when  time  is  zero,  and  thus,  the  initial 
values  of  the  nonlinear  ordinary  differential  equations  cannot 
be  computed. 

B.  RESULTS  FOR  THE  FLAT  PLATE  AND  THE  SPHERE 

Integral  equations  (1.39)  and  (1.43)  are  solved  using  the 
Runge-Kutta  method  of  orders  1,  3,  and  5.  The  first  order 
approximation  can  be  found  in  [5],  whereas  the  third  and  the 
fifth  order  approximations  are  described  in  sections  3(C)  and 
3(D),  respectively.  Solutions  of  the  nonlinear  ordinary 
differential  eq[uatlons  corresponding  to  (1.39)  and  (1.43)  are 
obtained  using  the  fourth-order  Runge-Kutta  method  developed 
by  Zurmuhl  [15].  The  results  show  that  solutions  of  a  high 
order  approximation  fall  below  those  of  a  lower  order 
approximation  (Fig.  5.1,  Fig.  5.2,  Fig  5.3). 
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rig.  S.l  Surface  temperature  of  a  Flat  Plate 
cooled  by  convection  and  radiation 
(Runce-Kutta  Method) . 


Fig.  yz  Surface  temperatture  of  a  Flat  Plate 
cooled  by  convection  and  radiation 
(Rtiace-Ctina  Method) . 
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Fig.  3.3  Surface  temperature  of  a  tphere 
cooled  by  convection  and  i^ation 
(Runge-Cutta  Method) . 
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with  respect  to  time  step,  we  do  not  have  the  sane  phenomenon 
as  in  the  order  of  approximation.  In  a  fixed  order 
approximation  method,  the  solution  curves  for  a  smaller  time 
step  fall  below  those  for  a  larger  time  step  at  small  times 
(approx,  less  than  0.2)  and  above  at  large  times  (Fig.  5.4). 


Eig.  3.4  Surface  temperature  of  a  flat  Plate 
cooled  by  convection  and  radiation 
(Runge-lCtitta  Method  of  the  First  Order) 

According  to  numerical  experiments,  the  stability  requirements 
for  Runge-Kutta  method  of  orders  1,  3,  and  5  are  approximately 
0.1,  0.01,  and  0.001,  respectively.  As  observed  earlier,  a 
drawback  of  the  Runge-Kutta  method  is  that  it  requires  the 
solution  of  a  heuristic  nonlinear  ordinary  differential 
equation  for  a  high  order  approximation.  This  leads  to  an 
attempt  to  use  an  easier  algorithm,  and  for  this  reason,  the 
finite  difference  method  was  implemented. 

Equations  (1.1)  and  (1.2a,b,c)  are  solved  for  the  flat 
plate  and  for  the  sphere.  The  extrapolation  formula  used  to 
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improve  the  accuracy  of  the  solutions  is  the  Aitken-Neville 
algorithm  (4.32).  The  results  for  various  time  steps  are 
presented  in  tables  1  and  2,  and  some  of  these  results  are 
plotted  in  Figures  5.5  and  5.6.  As  tables  1  and  2  show,  the 
situation  where  solution  curves  for  a  smaller  time  step  fall 
below  those  for  a  larger  time  step  holds  in  the  finite 


difference  method. 


table  1 


The  Finite  PifTerenee  Method  for  Various  Time  Steps 


Time 

At 

icr2 

ir® 

IQT^ 

0.01 

0.849395 

0.843059 

0.842539 

0.02 

0.797214 

0.793829 

0.793609 

0.764701 

0.762365 

0.762256 

0.04 

0.740419 

0.738609 

0.738559 

0.05 

0.720827 

0.719340 

0.719328 

0.06 

0.704316 

0.703048 

0.703063 

0.07 

0.689998 

0.688893 

0.688927 

0.08 

0.677333 

0.676351 

0.6764(X) 

0.09 

0.659595 

0.665075 

0.665137 

0.10 

0.655626 

0.654821 

0.654893 

0.20 

0.583564 

0583127 

0.30 

0.535891 

0.535564 

0.40 

0.496573 

0.496292 

0.50 

0.461310 

0.461056 

0.60 

0.428810 

0.428575 

0.70 

0.398614 

0.398394 

0.80 

0.370501 

0.370296 

0.344321 

0.344130 

liB 

0.319948 

0.319770 

The  Surfaea  at  a  Hat  llata  eeoled  by 

Cenvaetlaii  and  ftadiatioti. 
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tig.  5.S  Surface  temperature  of  a  flat  plate 
eooled  by  eonvaetion  and  radiation 
(finite  bifferenee  Method.  bUo.O  t) 


tabtl2 

The  finite  bifference  Method  far  Verlout  Time  Stept 


Time 

-At  1 

ir^ 

itr^ 

0.0i 

0.916249 

0.912694 

0.911001 

0.02 

0.882656 

0.880526 

0.860270 

0.858708 

0.842779 

0.841514 

0.05 

0.828114 

0.827055 

0.06 

0.814411 

0.07 

0.801988 

0.801114 

0.08 

0.791626 

0.792846 

0.09 

0.774505 

0.10 

0.775178 

0.20 

0.10 

0.656561 

0.40 

0.616906 

0.50 

0.584412 

0.60 

0557215 

0.70 

0.514052 

0.80 

0.514011 

0.90 

0.496514 

1.00 

0.481017 

The  surface  Temperature  of  t  iphere  cooled  by 
ConTectien  end  Radiation. 
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Fig.  S-6S\irfaee  temperattire  of  a  sphere 
cooled  by  convection  and  radiation 
(Finite  Difference  Method). 

Even  though  the  implicit  scheme  is  implemented  on  the 
boundary,  numerical  experiments  show  that  solutions  still 
exhibit  oscillation  when  a  large  time  step  was  chosen  (time 
step  >  0.01).  This  constraint  of  time  step  leads  to  large 
computational  times  for  large  time  solutions. 

Figures  5.7  and  5.8  show  representative  results  for  the 
Runge-Kutta  and  finite  difference  methods  where  a  flat  plate 
and  sphere  are  cooling.  Tables  3  and  4  show  that,  when  At  - 
0.01,  the  results  obtained  by  using  the  Runge-Kutta  method  of 
orders  2  and  3  compared  favourably  with  those  using  the  finite 
difference  method.  The  difference  of  the  solutions  by  using 
the  two  methods  is  less  than  3.1%  (relative  error)  in  average 
for  each  case. 
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Comparison  of  the  Runge-Kutta  and  the 
Finite  Difference  Methods  (At=0.01) 


Time 

1st  Order 

5rd  Order 

Finite 

DIff. 

0.10 

0.712925 

0.667017 

0.655626 

0.20 

0.619477 

0567586 

0.583564 

0.30 

0.562454 

0514915 

0.535891 

0.40 

0.516209 

0.477480 

0.496573 

0.50 

0.475091 

0.445843 

0.461310 

0.60 

0.437583 

0.417183 

0.428810 

0.70 

0.403124 

0.390595 

0.398614 

0.80 

0.371405 

0.365742 

0.370501 

0.90 

0.342189 

0.342459 

0.344321 

1.00 

0.315274 

0.320636 

0.319948 

The  Surface  Temperature  of  a  flat  Plate  cooled  by 
Convection  and  Radiation. 


Fic.  !1.7Comparrion  of  results  for  cooling 
a  flat  plate. 
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Tabte4 

Comparison  of  the  Runge-Kutta  and  the 
Finite  Difference  Methods  (At=0.01) 


Time 

1st  Order 

3fd  Order 

Finite 

Dlff. 

0.10 

0.785102 

0.754306 

0.775178 

020 

0.705908 

0.665482 

0.706513 

0.30 

0.651471 

0.622430 

0.656562 

0.40 

0.609637 

0591087 

0.616906 

0.50 

0.576138 

0.565188 

0.584412 

0.60 

0.548557 

0542988 

0.557215 

0.70 

0.525344 

0.523640 

0.534052 

0.80 

0.505453 

0.506582 

0.514033 

0.90 

0.488156 

0.491399 

0.496514 

1.00 

0.472926 

0.477772 

0.481017 

The  Surface  Temperature  of  a  Sphera  cooled  by 
Convection  and  Radiation. 


Fit.  3.8  Coffiparsion  of  rasultt  fbr  cooling 
of  a  sphera. 


79 


Lastly,  in  Figures  5.9  and  5.10,  the  results  of  two  special 
cases  solved  by  finite  difference  method  are  compared.  The 
graph  shows  that,  when  At  «  0.01,  the  surface  temperature  of 
a  flat  plate  fell  much  faster  than  that  of  a  sphere.  We 
believe  that  it  is  due  to  the  effect  of  the  boundary  condition 
at  X  =  0  and  the  difference  in  the  coefficient  of  the 
convective  term. 


Surface  Temp 
0.3 
0.25 
0.2 
0.15 
0. 1 
0.05 

rig.  5.9  Surface  temperature  of  a  flat  plate 
cooled  by  convection  and  radiation 
(Finite  Pifferenea  Method,  at =0.0 1) 


Fig.  5.10  Surface  temperature  of  a  sphere 
cooled  by  convection  end  radiation 
(Finite  Difference  Method,  at =0.0 1) 
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VI.  COMCLnSZOMS 


The  study  of  the  one  dimensional  heat  equation  subject  to 
combined  convective  and  radiative  boundary  conditions  in 
rectangular  coordinates  is  motivated  by  the  advent  of  space 
technology  where  ]cnowledge  of  the  temperature  of  bodies  in 
deep  space  is  necessary,  for  instance,  in  the  design  of  space 
shuttles . 

The  solids  are  assumed  to  be  homogeneous,  isotropic,  and 
opaque  to  thermal  radiation  and  to  have  temperature 
independent  physical  properties.  This  assumption  leads  to  a 
linear  heat  equation.  The  difficulty  of  the  problem  is 
determined  by  the  conditions  prescribed  at  the  boundaries. 
According  to  the  laws  of  physics,  the  heat  flux  of  the 
radiative  heat  transfer  is  proportional  to  the  fourth  power  of 
the  temperature  which  causes  nonlinearity  at  the  boundaries. 

Problems  of  this  type  are  first  solved  by  analytic 
techniques,  one  of  which  is  the  integral  transform  method.  In 
particular,  Laplace  transform  and  eigenvalue  expansion  are 
used.  The  solutions  which  are  explicitly  determined  at  the 
surface  for  two  special  cases,  neunely,  the  flat  plate  and  the 
sphere,  are  singular  nonlinear  Vol terra  integral  equations  of 
the  second  kind.  Although  they  are  not  practical  in 
determining  the  temperature  at  a  particular  time,  these 
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integral  equations  can  help  us  to  deduce  some  useful 
information  about  the  behaviour  of  the  surface  temperature. 

Since  the  analytic  solutions  found  for  the  problem  are  not 
practical  to  use,  numerical  techniques  are  considered  as  an 
alternative.  Two  numerical  schemes  that  are  used  to  deal  with 
the  resulting  integral  equation  are  the  Runge>Kutta  method  and 
the  successive  approximations  method.  Both  techniques  are 
studied  in  great  detail.  Ntunerical  solutions  show  that  the 
successive  approximations  method  is  "exact"  in  the  sense  that 
any  desired  accuracy  may  be  obtained  [3,4].  Additionally,  the 
closer  the  Initial  approximation  was  to  the  exact  solution, 
the  faster  the  method  of  successive  approximation  converged  to 
the  exact  solution.  Conditions  for  the  numerical  solution  and 
limitations  of  these  schemes  are  also  discussed. 

Another  numerical  technique  which  is  directly  applied  to 
the  governing  equations  is  presented  as  a  possible  alternative 
to  the  numerical  methods  previously  discussed.  It  is  the  well 
known  finite  difference  method  in  which  the  Crank-  Nicholson 
scheme,  the  backward  implicit  scheme,  and  the  Newton-Raphson 
method  are  combined  to  solve  for  the  surface  temperature. 

The  Runge-Kutta  methods  of  orders  1,  3,  and  5  are 
programmed  for  (1.39)  and  (1.43)  which  are  the  integral 
equations  corresponding  to  the  flat  plate  and  the  sphere, 
respectively.  The  numerical  results  are  presented  with 
respect  to  their  orders  and  to  their  time  steps.  The  data 
reveal  the  following  phenomena.  First,  the  solutions  of  a 
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high  order  approximation  fall  below  those  of  a  lower  order 
approximation.  This  phenomena  is  a  result  of  the  higher  order 
approximations  closing  in  on  the  actual  solution.  Second,  the 
first  phenomenon  does  not  occur  in  the  solutions  for  various 
time  steps  with  a  fixed  approximation  order.  The  main  result 
here  is  that  a  smaller  step  size  determines  the  surface 
temperature  for  very  small  times  (0£tj:0.2)  more  accurately  and 
a  larger  step  size  determines  the  surface  temperature  for 
larger  times  (tzO.2)  more  accurately.  Third,  the  agreement 
between  the  1—,  3—,  and  5—  order  Runge-Kutta  approximations 
is  better  for  the  sphere  than  that  for  the  plate.  Physically, 
this  is  due  to  the  fact  that  the  boundary  surface  area  to 
total  volume  ratio  is  largest  for  the  sphere  and  smallest  for 
the  plate.  The  reason  for  this  trend  is  that  the  larger  the 
ratio  the  more  uniform  will  the  temperature  be  throughout  the 
body.  The  fourth  phenomenon  is  that  the  accuracy  of  the 
approximation  increases  with  time.  For  large  values  of  time, 
the  rate  of  change  of  temperature  is  reduced,  as  would  be 
expected  from  the  influence  of  the  fourth  power  term  (U^)  . 
Since  the  Runge-Kutta  method  did  not  offer  any  efficiency  in 
the  area  of  high  order  approximations,  the  finite  difference 
method  was  considered. 

Equations  (1.1)  and  (1.2a,b,c}  are  solved  numerically 
using  the  finite  difference  method  for  both  the  flat  plate  and 
the  sphere.  The  results  for  various  time  steps  are  presented. 
The  table  shows  that  the  second  phenomenon  found  in  the  Runge- 
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Kutta  method  again  occurs  in  the  solutions  generated  by  the 
finite  difference  method  with  respect  to  time  step.  Again,  as 
in  the  Runge-Kutta  method  the  smaller  step  size  determines  the 
surface  temperature  more  accurately  for  small  times  and  the 
larger  step  size  determines  the  surface  temperature  more 
accurately  for  large  times. 

Finally,  two  comparisons  are  made  of  the  numerical 
solutions.  The  first  is  of  the  Runge-Kutta  method  and  the 
finite  difference  method.  The  results  show  that  there  is  a 
good  agreement  between  the  two  methods,  and  the  difference 
between  their  solutions  are,  on  the  average,  less  than  3.1%  in 
both  cases.  The  second  comparison  was  made  between  the 
solution  of  a  flat  plate  and  that  of  a  sphere.  The  finite 
difference  method  conveys  that  temperature  of  a  flat  plate 
decays  much  faster  than  that  of  a  sphere.  This  result  was 
expected  for  the  transient  heat  conduction  with  linear 
boundary  conditions.  This  could  be  due  to  a  larger  area  on 
the  plate  exposed  to  the  uniform  boundary  layer. 
Additionally,  this  result  could  be  caused  by  the  sphere  having 
a  larger  surface  area  to  volume  ratio;  thus  the  sphere  would 
have  a  more  uniform  temperature  distribution  throughout  the 
body  resulting  in  a  slower  decay  of  surface  temperature. 

Comments  of  a  more  general  nature  are  included. 

1.  The  convection  mode  of  heat  transfer  appears  to  be 
dominant  as  the  dimensionless  temperature  approaches 
uniformity  for  a  plate  cooling  to  a  zero  environment.  This 
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result  is  due  to  the  fact  that  U*  is  approaching  zero  at  much 
a  faster  rate  then  U. 

2.  Physically,  the  adiabatic  or  initial  temperature 
cannot  be  equal  to  absolute  zero,  however,  in  many  situations 
the  temperature  ratio  of  adiabatic  surface  temperature  to 
initial  temperature  can  be  very  small. 

3.  For  cooling  and  heating  the  solutions  are  initially 
inaccurate  due  to  the  fact  that  at  t=0  the  linearized  heat 
flux  is  not  equal  to  the  actual  flux. 

4 .  For  a  set  time  step  size  the  number  of  iterations 
required  to  meet  a  set  accuracy  is  determined  by  which 
surface  is  receiving  the  highest  heat  rate. 

5.  The  time  required  to  achieve  a  particular  surface 
temperature  during  cooling  decreases  as  the  ratio  of  the 
environment  temperature  to  initial  solid  temperature 
increases . 

To  conclude  this  thesis,  a  numerical  scheme  is  proposed  as 
an  alternative  to  the  existing  ntuaerical  methods.  The  method 
of  successive  approximations  is  described  in  Chapter  II.  One 
of  the  major  difficulties  of  that  method  is  choosing  the 
initial  approximation  for  the  iteration  procedure.  As 
mentioned  earlier,  the  convergence  of  the  algorithm  can  be 
accelerated  if  one  could  obtain  an  initial  approximation  which 
is  close  to  the  exact  solution.  To  determine  this  value,  one 
could  first  use  the  finite  difference  method  (without  the 
extrapolation  algorithm)  described  in  Chapter  IV  to  determine 
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the  surface  temperature.  Then,  by  treating  it  as  an  initial 
approximation,  the  method  of  successive  approximations  is 
applied  to  obtain  the  solution.  We  believe  that  the 
temperature  obtained  by  using  the  finite  difference 
approximation  for  the  exact  solution  would  be  a  better 
solution  than  the  temperature  at  the  previous  time  level.  In 
addition,  this  technique  would  allow  larger  time  steps. 
However,  as  far  as  the  author  is  aware,  nothing  has  been 
proved  for  this  method,  and  the  analytical  and  numerical 
justifications  for  the  algorithm  are  left  open. 
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APPBMDIX-A 


To  provide  a  better  understanding  of  the  results  for  Ii(t)  and 
Ik(t)  corresponding  to  (3.29)  and  (3.30)  respectively, 
consider  the  following  exeunple  where  N  3  in  (3.17): 


Pi*i  Pa*!  Pa*! 

Iiit) 

-Pi*!  -Pa*!  -p3*! 

liit) 

= 

Ba 

Pi*!  Pa*!  Pa*! , 

JaCt) 

^3 

where 

Px«!  P2«l  p3«r 

•A  =  -Pi«i  -Pa«2  ~P3«3  •  (A.  2) 

Pi«i  Pa*!  Pa*! . 

Using  Cramer's  rule 


(A. 3) 


where 
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Pl«!  P2«i  P3«l 

Dec (A)  =  -pittl  -P2O2  -Patt?  . 

Pi4  Pjal  Pa®* 


(A. 4) 


Using  a  fundamental  property  of  determinants  (A. 4)  can  be 
written  as 


111 


Det{A)  =  -P1P2P3 (aia2a2)^i  • 


*  e,* 

1  «2  *3 


(A. 5) 


Column  reduction  gives 


Dec  (A)  =  -PiP2P3(a?aia^)  «i 

«!  (aj-ai)  («3-ai) 


(A. 6) 


Define 


a2  -(aa-tti)  . 

«!  (a2-«i)  (a|-a?) 


(A. 7) 


then 
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ret  (A') 


-(a|-a?)  -ial-al) 
(a^-aj)  (ffs-aj) 


(A.  8) 


By  a  similar  manipulation  (A. 3)  becomes 


Ji(t) 


Pitti  ret{A0 


(A. 9) 


After  working  through  a  considerable  amount  of  algebra  both 
(A. 3)  and  (A. 9)  give 

j  (t)  =  *  B^[a\  *  al)  * 

^  ^xallalal  -  alial  +  a|)  +  aj] 

In  a  similar  manner  IjCt)  and  IjCt)  can  be  determined.  Without 
loss  of  generality  Ik(t)  can  also  be  found. 
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APPEHDIZ-B 

Th«  TlioBas  Algoritlm 


The  equations  are: 

where  1  ^  1  ^  N  with  ax  =  =  0,  and  N  is  the  number  of  nodes 

in  the  domain. 

The  algorithm  is  as  follows: 

DO  10  i  *  2,N 
ratio  =  ax/bi 

bi  =  bi  -  ratio  *  Ci.x 
di  =  di  -  ratio  *  dx-i 

10  CONTINUE 

Un  *  dn/bN 

DO  20  i  =  N-1,1 

Ui  =  (dj  —  Cl  *  di+i)/bi 

20  CONTINUE 
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Your  participation  in  this  research  and  response  by  12 
March  1993  is  greatly  appreciated.  If  you  need  cuiy  additional 
details,  contact  LCDR  Richard  Mendez  (408)  759-9783/LT  Gerald 
Rivas  (408)  655-1625,  or  by  writing  to: 


LT  Gerald  A.  Rivas,  SC,  USN 
SMC  #2715 

Naval  Postgraduate  School 
Monterey,  CA  93943-5000 


Name: 
Address : 
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ALLOTOOmT 


An  authorization  by  the  head  (or  other  authorized  employee)  of 
an  operating  agency  which  assigns  a  specified  amount  of  money 
to  subordinate  units.  The  amount  allotted  by  the  agency 
cannot  exceed  the  amount  apportioned  by  the  Office  of 
Management  and  Budget  (0MB) . 

Synonyms :  None 
Antonyms :  None 


Do  you  agree  with  this  definition? 


. 1 . 2 . 3 . . 4-- . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 


COMMENT: 


Synonyms  t 
Antonyms: 
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AWARD 


(1)  The  formal  acceptance  of  an  offeror's  bid  or  proposal. 

(2)  Notification  of  intent  to  give  a  contract. 

(3)  Transmittal  of  advance  authorization  to  proceed  (e.g. 
letter  contract) . 

Synonyms :  None 
Antonyms :  None 


Do  you  agree  with  this  definition? 

. 1 . 2 . -3 . 4 . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 


COMMENT: 


Synonyms : 
Antonyms  t 
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The  formal  act  of  acknowledging  in  writing  and  affirming  by 
signature  that: 

-  some  act  has  or  has  not  been  performed; 

-  some  event  has  or  has  not  occurred; 

-  some  legal  formality  has  or  has  not  been  complied  with;  or 

-  some  condition  exists  or  does  not  exist. 

Synonyms  t  None 
Antonyms :  None 


Do  you  agree  with  this  definition? 


. 1 . 2 . 3 . 4 . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 


COMMENT: 


COWSULTAMT 

A  person  having  specialized  education  and/or  broad  experience 
which  uniquely  qualifies  them  to  be  called  upon  to  furnish 
expert  advice  on  highly  specialized  matters  and  recommend 
solutions  to  particular  problems. 

Synonyms:  Advisor,  Expert,  Subject  Matter  Expert  (SME) 
Antonyms  t  None 


Do  you  agree  with  this  definition? 


. 1 . 2 . 3- . 4 . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 


COMMENT: 


Synonyms : 
Antonyms : 
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COST  OBJECTIVE 

(1)  A  function,  organizational  subdivision,  contract,  or 
other  work  unit  for  which  cost  data  is  desired  and  for  which 
provision  is  made  to  accumulate  and  measure  the  cost  of 
processes,  products,  capitalized  projects,  and  so  forth. 

(2)  Cost  goal  established  for  the  completion  of  an  element 
of  work. 

(3)  Goal  established  for  contract  cost  to  be  achieved  during 
contract  negotiations. 

Synonyms:  Cost  Center,  Cost  Goal,  Target  Cost 

Antonyms :  None 


Do  you  agree  with  this  definition? 


. 1 . 2 . 3 . . 4 . 5--- 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE 
AGREE  RESERVATION  RESERVATION 


- 6 - 

STRONGLY 

DISAGREE 


COMMENT; 


Synonyms : 
Antonyms : 


09 


BSCALATION 


(1)  A  term  traditionally  used  to  indicate  an  upward  or  (more 
rarely)  a  downward  movement  of  price.  "Economic  Price 
Adjustment"  is  the  contemporary  term  used  to  express  the 
application  of  escalation  by  specified  procedures. 

(2)  In  Government  contracting  refers  to  an  amount  or  percent 
by  which  a  contract  price  may  be  adjusted  if  predefixied 
contingencies  occur,  such  as  changes  in  the  vendor's  raw 
material  costs  or  labor  costs.  The  amount  of  the  "escalation" 
is  usually  tied  to  some  predetermined  price  index. 

Synonyms:  Economic  Price  Adjustment 
Antonyms :  None 


Do  you  agree  with  this  definition? 


. 1 . 2 . 3 . ---4 . 5 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE 
AGREE  RESERVATION  RESERVATION 


- 6 - 

STRONGLY 

DISAGREE 


COMMENT: 


Synonyms : 
Antonyms : 
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BXPBWSB 


Costs  of  operation  and  maintenance  of  activities  on  the 
accrual  basis  for  a  fiscal  period,  as  distinguished  from 
capital  costs  that  will  be  depreciated  over  their  approximate 
service  life. 

Synonyms  t  Costs 
Antonyms :  Revenue ,  Income 


Do  you  agree  with  this  definition? 


- 1 . 2 . 3 . -4--- . 5-.- 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE 
AGREE  RESERVATION  RESERVATION 


- 6 - 

STRONGLY 

DISAGREE 


COMMENT: 


Synonyms : 
Antonyms : 
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IKDUSTRIAL  PLANT  BQUlPKgMT  (IPg) 

Plant  equipment  acquired  by  the  Government,  exceeding  an 
established  acquisition  cost  threshold,  used  for  the  purpose 
of  cutting,  edjrading,  grinding,  shaping,  forming,  joining, 
testing,  measuring,  heating,  treating  or  otherwise  altering 
the  physical,  electrical  or  chemical  properties  of  materials, 
components  or  other  end  items  entailed  in  manufacturing, 
maintenance,  supply,  processing,  assembly  or  research  and 
development  operations . 

Synonyms :  None 
Antonyms :  None 


Do  you  agree  with  this  definition? 


. 1 . 2 . 3 . 4 - - 5 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE 
AGREE  RESERVATION  RESERVATION 


- 6 - 

STRONGLY 

DISAGREE 


COMMENT: 


Synonyms : 
Antonyms  t 
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QEEJSEXS 


A  cost  balancing  action  whereby  a  claim  may  be  canceled  or 
lessened  by  a  counterclaim. 

Defective  pricing:  Allowable  understatements  (e.g., 
counterclaims  or  cost  proposal  errors  that  are  favorable  to 
the  contractor)  which  are  reduced  by  overstatements  of  cost 
that  arise  under  a  defective  pricing  case.  In  order  to 
eliminate  an  increase  in  the  contract  price  the  offset  cannot 
exceed  the  extent  of  the  overstatement. 

Administrative  Offset:  A  procedure  to  collect  a  debt  owed 
to  the  Government  by  withholding  money  payable  to  contractor 
under  a  contract,  in  order  to  satisfy  the  contractor's  debt 
which  arose  independently  of  that  contract  and  which  are  in 
compliance  with  the  Federal  Claims  Collection  Act  of  1966. 

Synonyias:  Counterclaim,  Setoff 
Antonyms :  None 


Do  you  agree  with  this  definition? 


. 1 . 2 . 3 . 4 . 5 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE 
AGREE  RESERVATION  RESERVATION 


- 6 - 

STRONGLY 

DISAGREE 


COMMENT: 


Synonyms  t 
Antonyms : 
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WBIGBTHD  AVBRAGB  COST  METHOD 

A  method  cf.  determining  the  average  unit  cost  of  inventory  and 
by  implication  an  aid  in  determining  the  cost  of  goods  made, 
sold,  or  held  for  future  sale  or  incorporation  into  higher 
level  end  items.  Under  this  technique,  costs  are  periodically 
computed  by  adding  the  sum  of  the  costs  of  beginning  inventory 
with  the  sum  of  the  costs  of  subsequent  purchases  and  dividing 
by  the  total  number  of  units. 

Synonyms  t  None 
Antonyms :  None 


Do  you  agree  with  this  definition? 


. 1 . 2 . 3 . 4 - - 5  - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE 
AGREE  RESERVATION  RESERVATION 


- 6 - 

STRONGLY 

DISAGREE 


COMMENT: 


Synonyms : 
Antonyms : 
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Allocation 


Funding:  An  amount  of  money  transferred  from  one  agency, 
bureau  or  account  that  is  set  aside  in  an  appropriation  of  the 
various  committees  having  spending  responsibilities  to  carry 
out  the  purposes  of  the  parent  appropriation  or  fund. 

Within  DOD,  the  money  is  being  transferred  from  the  services 
to  the  appropriate  MAJCOMS. 

Financial:  A  cost  accounting  procedure  which  results  in  a 
reasonable  distribution  of  costs  cunong  one  or  more  cost 
cbjectives  (e.g.,  products,  programs,  contracts,  and 
activities)  .  This  includes  both  direct  assignment  of  costs  and 
the  reassignment  of  a  share  from  an  indirect  pool. 

Synonyms :  None 
Antonyms  t  None 


Do  you  agree  with  this  definition? 


. 1 . 2 . 3 . 4 . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 


COMMENT: 


Synonyms : 
Antonyms : 
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j&wQrttgfttir0.p 

The  systematic  reduction  of  an  indebtedness  or  recorded  asset 
value  over  a  specific  period  of  time  by  periodic  payments  to 
a  creditor  or  charges  to  an  expense,  in  accordance  with 
generally  accepted  accounting  procedures  or  principles. 

Synonyms t  Liquidation,  Allocation,  Writeoff 
Antonyms:  Direct  charge 


Do  you  agree  with  this  definition? 


. 1 . 2 . 3 - 4 . 5 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE 
AGREE  RESERVATION  RESERVATION 


- 6 - 

STRONGLY 

DISAGREE 


COMMENT: 


Synonyms : 
Antonyms: 
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Bid 

An  offer  to  perform  a  contract  by  providing  lcd)or  and  or 
material  for  a  specific  price.  In  federal  government 
contracting,  this  offer  is  provided  in  response  to  an 
invitation  for  bid. 


Synonyms  t  None 
Antonyms  t  None 


Do  you  agree  with  this  definition? 


. 1 . 2 . 3 . . -4 . . 5--- . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 


COMMENT: 


Synonyms : 
Antonyms : 
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The  act  by  an  authorized  individual  affirming  the  intent  of  an 
agency  or  company  to  take  or  accept  a  defined  action  not  yet 
formalized  by  execution  of  a  contract. 

Funding:  A  firm  administrative  reservation  of  funds  based  upon 
firm  procurement  directions,  orders,  requisitions,  certified 
purchase  requests,  and  budgetary  authorizations  which  set 
aside  certain  funds  for  a  particular  contract  without  further 
recourse  to  the  official  responsible  for  certifying  the 
availcibility  of  funds. 

Within  DOD,  reservation  of  funds  are  set  aside  by  the 
appropriate  operating  division  (wing  or  base)  for  use  on  a 
particular  item. 

Accounting:  The  method  of  accounting  for  the  available  balance 
of  an  appropriation,  fund,  or  contract  authorization  whereby 
commitments  are  recorded  in  the  accounts  as  reductions  of  the 
available  balance. 

Synonyms :  None 
Antonyms :  None 

Do  you  agree  with  this  definition? 


. 1 . 2 . 3 . 4 . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 


COMMENT: 


Synonyms  t 
Antonyms  t 
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Cost 

(1)  For  the  Seller:  The  amount  of  money  or  equivalent 

incurred  for  supplies  or  services  exclusive  of  profit  or  fee. 

(2)  For  the  Buyer:  The  amount  of  money  or  equivalent  paid 
for  supplies  or  services  including  profit  or  fee. 

Synonyms:  Expense,  Consideration,  Charge,  Total  Cost 

Antonyms :  None 


Do  you  agree  with  this  definition? 


. 1 . 2 . 3 . 4 . 5 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE 
AGREE  RESERVATION  RESERVATION 


- 6 - 

STRONGLY 

DISAGREE 


COMMENT: 


Synonyms : 
Antonyms : 
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Dalipqp^gn<?y 

(1)  Failure,  omission,  or  violation  of  contractual 
obligation  or  duty. 

(2)  The  actual  failure  by  the  contractor  to  meet  the 
contract  delivery  or  performance  schedule,  or  the  potential 
failure  to  do  so  by  failing  to  maintain  required  progress  in 
contract  performance  as  required  by  the  contract  delivery  or 
performance  schedule 


Synonyms:  Overdue,  Tardy,  Late 
Antonyms:  Early,  Accelerated,  Timely 


Do  you  agree  with  this  definition? 


. 1 . 2 . 3--- . 4 . 5 .  -6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 

COMMENT; 


Synonyms : 
Antonyms : 
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Executed  Contract 

A  written  document  which  has  been  signed  by  both  parties  and 
itiailed  or  otherwise  furnished  to  each  party,  which  expresses 
the  requirements,  terms,  and  conditions  to  be  met  by  each 
party . 


Synonyms :  None 
Antonyms :  None 


Do  you  agree  with  this  definition? 


. 1 . 2 . 3- . -4 . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 


COMMENT: 


Synonyms : 
Antonyms : 
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Indwwdeat  Cogt 

A  cost  estimate  developed  outside  the  normal  advocacy 
channels,  independent  of  any  cost  information  provided  by  the 
offeror,  used  for  the  purpose  of  cort^aring  with  bids  or 
proposals.  Preparation  of  independent  costs  estimates 
generally  include  representations  from  the  areas  of  cost 
analysis,  procurement,  production  management ,  engineering,  and 
program  management. 

Synonyms:  Independent  Government  Cost  Estimate  (IGCE) 
Antonyms :  None 


Do  you  agree  with  this  definition? 


. 1 . 2 . 3 . ---4 . . 5--- 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE 
AGREE  RESERVATION  RESERVATION 


- 6 - 

STRONGLY 

DISAGREE 


COMMENT: 


Synonyms : 
Antonyms : 
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Novation  Agreement 

A  legal  instrument,  executed  by  the  parties  to  a  contract  and 
a  successor  in  interest,  which  transfers  all  obligations  and 
rights  under  the  contract  to  the  successor. 

The  government  may  recognize  a  third  party  as  a  successor  of 
a  government  contract  when  the  third  party's  interests  arises 
out  of  the  transfer  of  1)  all  the  contractor's  assets,  or  2) 
the  entire  portion  of  the  assets  involved  in  the  performing  a 
contract . 

Synonyms :  None 
Antonyms :  None 


Do  you  agree  with  this  definition? 


. 1 . 2 . 3 . 4 . 5-  - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE 
AGREE  RESERVATION  RESERVATION 


- 6 - 

STRONGLY 

DISAGREE 


COMMENT: 


Synonyms : 
Antonyms : 
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BpYAlty 

Con5)ensation  paid  to  the  owner,  vendor  or  lessor  of  personal, 
real,  tangible  or  intangible  property  for  the  use  of  that 
property.  Usually  a  percentage  of  the  selling  price  of  goods 
and  services,  production  of  which  employs  the  property. 

Synonyms:  Commission  Payment,  Use  Fee 
Antonyms:  Royalty  Free  Use 


Do  you  agree  with  this  definition? 


. 1 . 2 . 3 - 4 . ---5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 


COMMENT; 


Synonyms : 
Antonyms : 
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APPENDIX  C:  FOLLON-UP  QUESTIONNAIRE 


This  appendix  represents  the  follow-on  questionnaire  which 


was  mailed  to  the  respondents  from  the  initial  questionnaire 


who  provided  neunes  and  addresses. 

A.  FOLLOW- ON  QUESTIONNAIRE 

Thank  you  for  coir^leting  the  initial  questionnaire  and 
participating  in  this  follow-on  research  to  arrive  at 
consensus  definitions  of  contracting  terminology.  Your 
efforts  have  provided  an  excellent  base  for  the  establishment 
of  a  consensus.  This  questionnaire  will  only  be  sent  to  those 
who  responded  to  the  Initial  questionnaire,  so  your  continued 
participation  is  very  inportant. 

As  a  reminder:  Graduate  students  at  the  Naval 
Postgraduate  School,  Monterey,  California,  and  the  Air  Force 
Institute  of  Technology,  Wright -Patterson  Air  Force  Base, 
Ohio,  are  conducting  research  to  derive  baseline  definitions 
for  commonly  used  acquisition  words  or  phrases.  When  the 
project  is  conplete,  the  definitions  will  be  included  as  part 
of  a  professional  dictionary  of  contracting  terminology  that 
will  be  published  by  the  NCMA.  The  purpose  of  the  dictionary 
is  two  fold.  First,  to  provide  an  educational  tool  to  those 
unfamiliar  with  the  acquisition  process.  Second  to  provide  a 
reference  document  for  those  working  in  the  field.  This 
research  is  an  ongoing  effort  in  obtaining  feedback  from 
contracting  professionals  regarding  proposed  definitions  of 
contracting  terms.  It  differs  from  the  previous  research  in 
that  it  is  taking  terms  from  previous  efforts  which  generated 
significant  diversity,  and  is  refining  them  using  the  Delphi 
Technique.  All  terms  were  synthesized  from  collected 
definitions.  Government  regulations  and  contracting  literature 
and  were  reviewed  once  by  NCMA  Fellows  and  Certified 
professionals  prior  to  your  input  on  the  initial 
questionnaire . 

Attached  for  your  review  are  the  revised  definitions  and 
selected  comments  from  the  initial  questionaires .  The 
definitions  were  revised  by  the  researchers  and  reviewed  by  a 
committee  of  contracting  professionals  for  coitpliance  with  the 
consensus.  Please  review  the  revised  definitions  and  indicate 
your  agreement  level  on  the  scale  provided  from  1  to  6.  If 
you  have  any  disagreements  or  comments,  please  either  annotate 
them  where  applicable,  or  write  them  on  the  space  provided. 
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Your  continued  participation  in  this  research  and  response  by 
26  April  1993  is  greatly  appreciated.  If  you  need  any 
additional  details,  contact  LCDR  Richard  Mendez  or  LT  Gerald 
Rivas  by  telephone  at  (408)  656-2536  (Administrative  Science 
Curriculum  Office) ,  or  by  writing  to: 

LT  Gerald  A.  Rivas 
SMC  #2715 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 
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Original  Definition: 

ALLonrorr 


An  authorization  by  the  head  (or  other  authorized  employee)  of 
an  operating  agency  which  assigns  a  specified  amount  of  money 
to  subordinate  units.  The  amount  allotted  by  the  agency 
cannot  exceed  the  amount  apportioned  by  the  Office  of 
Management  and  Budget  (0MB) . 

Synonyms :  None 
Antonyms  t  None 
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Comments t 

Allotments  can  be  made  by  other  than  "heads"  to  "subordinate 
units" . 

Allotments  go  farther  than  to  subordinate  agency  units.  They 
end  up  being  made  to  programs /projects  and  individual 
contracts . 
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Periodicity  of  allotments,  i.e.  quarterly /annually. 

Add  to  end  of  first  sentence  ",  projects  or  activities." 

Can  negotiation  go  on  between  subordinate  units. 

Government  Contracts  Reference  Book  definition:  "In  OOD,  the 
process  by  which  commanders.  Major  Commanders,  or  Special 
Operating  Agencies  distribute  their  allocated  funds  to 
themselves,  to  installation  commanders  or  to  other  subordinate 
organizations.  This  process  may  continue  into  as  many  sub 
allotments  as  necessary." 

Synonyms:  Funding,  Budgeted  Amount,  Obligation, 
Appropriation,  Public  Troth. 

Antonyms : 

Revised  Definition: 

AT.T.fVrMgMT 

An  authorization  by  the  head  (or  other  authorized  employee)  of 
an  operating  agency  which  assigns  a  specified  amount  of  money 
to  subordinate  units,  projects  or  activities.  The  amount 
allotted  by  the  agency  cannot  exceed  the  amount  apportioned  by 
the  Office  of  Management  and  Budget  (0MB) . 

Synonyms :  Funding . 

Antonyms :  None 

Do  you  agree  with  this  definition? 

. 1 . 2 . 3 . 4 . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 

C0M1ENT ;  _ 


SYNONYMS: 

ANTONYMS: 
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Original  Definition 
AWARD 


(1)  The  formal  acceptance  of  an  offeror's  bid  or  proposal. 

(2)  Notification  of  intent  to  give  a  contract. 

(3)  Transmittal  of  advance  authorization  to  proceed  (e.g. 
letter  contract) . 

Synonyms :  None 
Antonyms  t  None 

Survey  Results 
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Comments : 

In  (2)  change  "give"  to  "let". 

Very  essential  to  establish  clear  and  concise  communication 
between  contractor  and  Government  contracting  personnel. 

Delete  (2)  .  Notification  of  intent  to  award  does  not 
constitute  award.  Enhance  definition  by  adding  "Notice  of 
Award,  Notice  to  Proceed"  to  (3) . 
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Award  is  the  process  through  which  a  buyer  and  seller  come  to 
agreement  over  the  terms  of  a  contract,  "award"  is  always 
rendered  by  the  buyer. 

Item  (2)  is  too  broad  -  needs  more  specifics  -  leaves  out 
FAR/DFAR. 

Item  (3)  is  a  conditional  award. 

Item  (2)  may  or  may  not  result  in  a  contract  depending  on 
state  of  negotiation. 

"Notification  of  intent"  and  "advance  authorization"  are  not 
considered  "award" . 

Item  (1)  needs  to  mention  "formal  acceptance  by  an  authorized 
official  of  the  Government". 

Item  (2) ,  intent  is  not  an  award. 

Items  (2)  and  (3)  are  a  result  of  Item  (l)  and  they  infer 
acceptance . 

Add  to  end  of  (1)  "as  offered." 

Item  (2)  requires  return  notice  of  acceptance  by  offeror. 

Reference  Book  Definition:  "The  notification  by  the 

Government  that  it  will  contract  with  a  private  party.  The 
award  of  a  contract  is  usually  made  by  Acceptance  of  an  Offer 
that  has  been  made  by  an  offeror.  In  procurements  by  sealed 
bidding,  the  contracting  officer  makes  a  contract  award  by 
written  notice,  within  the  time  for  acceptance  specified  in 
the  bid  or  extension,  to  the  responsible  bidder  whose  bid, 
conforming  to  solicitation,  is  the  most  advantageous  to  the 
Government,  considering  only  price  and  price- related  factors 
included  in  the  solicitation.  In  procurements  by  negotiation, 
the  contracting  officer  awards  a  contract  with  reasonedDle 
promptness  to  the  successful  offeror  (the  source  whose  best 
and  final  offer  (BAFO)  is  most  advantageous  to  the  Government 
considering  price  and  other  factors  included  in  the 
solicitation)  by  transmitting  a  written  notice  of  award  to 
that  offeror. 

Synonyoisx  Contract,  Win,  Def initization  of  Contract. 
Antonyms:  Loss 
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Revised  Definition 


AWARD 

(1)  The  formal  acceptance  of  an  offeror's  bid  or  proposal. 

(2)  Transmittal  of  advance  authorization  to  proceed  (e.g. 
letter  contract) . 

Synonyms  t  None 
Antonyms :  None 


Do  you  agree  with  this  definition? 

. 1 . 2 . 3 - - -4 . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 

COMMENT  : 


SYNONYMS: 

ANTONYMS; 
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Original  Definition! 


gBRTIfygATlOW 

The  formal  act  of  acknowledging  in  writing  and  affirming  by 
signature  that: 

■  some  act  has  or  has  not  been  performed; 

-  some  event  has  or  has  not  occurred; 

-  some  legal  formality  has  or  has  not  been  coii:5)lied  with;  or 

-  some  condition  exists  or  does  not  exist. 

Synonyms :  None 
Antonyms :  None 


Survey  Results 
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Comments  t 

Needs  to  follow  FAR/DFARS/CPR  more  closely. 

Need  to  mention  legal  accountability  of  certifier. 

Reword  definition  as  a  positive  statement.  Strike  "or  has 
not"  and  "or  does  not"  from  sentences. 


Requirement  should  be  mentioned  in  definition. 
Delete  item  3 . 


Synonyms : 

Antonyms : 

Revised  Definition: 

CBRTIPICATIQM 

The  formal  act  of  acknowledging  in  writing  and  affirming  by 
signature  that: 

-  some  act  has  or  has  not  been  performed; 

-  some  event  has  or  has  not  occurred; 

-  some  legal  formality  has  or  has  not  been  coital ied  with;  or 

-  some  condition  exists  or  does  not  exist. 

Synonyms :  None 
Antonyms :  None 

Do  you  agree  with  this  definition? 

. 1 . 2 . 3 . 4 . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 

COMMENT  ; 


SYNONYMS: 

ANTONYMS: 
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Original  Definition: 


1 


CQiTSTff.T|tffr 

A  person  having  specialired  education  and/or  broad  experience 
which  uniquely  qualifies  them  to  be  called  upon  to  furnish 
expert  advice  on  highly  specialized  matters  and  recommend 
solutions  to  particular  problems. 

Synonyms:  Advisor,  Expert,  Subject  Matter  Expert  (SME) 
Antonyms :  None 
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Survey  Results 

Consu 1 tant 

First  Round 


Comments : 

Change  "them"  to  "him/her". 

Change  "broad"  to  "extensive". 

Consultants  may  be  called  upon  to  furnish  guidance  or  to 
advocate,  for  the  benefit  of  an  individual  or  entity. 

Add  after  "expert  advice"  -  "or  opinions". 
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Add  after  "called  upon" 


"by  the  Federal  Government". 


Add  to  end  of  definition  "of  a  non- inherently  governmental 
nature . " 

Definition  too  narrow.  Is  deliverable  required? 

Consultant  connotes  a  business  relationship  unlike  synonyms. 

Change  "highly  specialized"  to  "various  or  relevant". 

Add  to  end  of  definition  "and  is  so  called  upon  for  that 
specific  purpose."  this  will  exclude  persons  already 
obligated  by  govt  contract. 

Synonyms s  Specialist,  Facilitator,  Authority. 

Antonyms:  Employee. 

Revised  Definition: 

CONSULTANT 

A  person  having  specialized  education  and/or  broad  experience 
which  uniquely  qualifies  him/her  to  be  called  upon  to  furnish 
expert  advice  or  opinions  on  highly  specialized  matters  and 
recommend  solutions  to  particular  problems. 

Synonyms:  Advisor,  Expert,  Svibject  Matter  Expert  (SME)  , 

Specialist,  Authority. 

An  tonyms :  None 

Do  you  agree  with  this  definition? 

. 1 . 2 . 3 . . 4 . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 

COMMENT  : 


SYNONYMS : 
ANTONYMS: 
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Original  Dafinitiont 


COST  OBJBCTIVB 

(1)  A  function,  organizational  subdivision,  contract,  or 
other  work  unit  for  which  cost  data  is  desired  and  for  which 
provision  is  made  to  accumulate  and  measure  the  cost  of 
processes,  products,  capitalized  projects,  and  so  forth. 

(2)  Cost  goal  established  for  the  completion  of  an  element 
of  work. 

(3)  Goal  established  for  contract  cost  to  be  achieved  during 
contract  negotiations. 

Synonyms:  Cost  Center,  Cost  Goal,  Target  Cost 
Antonyms :  None 


Survey  Results 
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Comments t 

Delete  (1)  and  "Cost  Center"  synonym. 

Add  to  (1)  after  "data  is  desired"  -  "and/or  required". 
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Change  definition  to  "Cost  objective  is  a  measure  of 
applicable  dollars  to  a  defined  task/work  effort.  Can  apply 
to  a  contract,  organization  or  other  work  unit." 

Synonyms:  Cost  Segment. 

Antonyms i  None 

Revised  Definition: 

COST  OBJBCTIVB 

Accounting:  A  function,  contract,  or  other  work  unit  for 

which  cost  data  is  desired  and  for  which  provision  is  made  to 
accumulate  and  measure  the  cost  of  processes,  products, 
capitalized  projects,  and  so  forth. 

Program  Management :  Cost  goal  established  for  the  completion 
of  an  element  of  work. 

Negotiations:  Goal  established  for  contract  cost  to  be 

achieved  during  contract  negotiations. 


Synonyms:  Cost  Goal,  Target  Cost 
Antonyms :  None 

Do  you  agree  with  this  definition? 


. 1 . 2--- 

STRONGLY  AGREE 

. 3 . 

AGREE  W/ 

. 4 . 

DISAGREE  W/ 

- 5 . 

DISAGREE 
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RESERVATION 

RESERVATION 

DISAGREE 
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SYNONYMS: 

ANTONYMS; 
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Original  Definition; 


KSCAIATIOW 

(1)  A  term  traditionally  used  to  indicate  an  upward  or  (more 
rarely)  a  downward  movement  of  price.  "Economic  Price 
Adjustment"  is  the  contemporary  term  used  to  express  the 
application  of  escalation  by  specified  procedures. 

(2)  In  Government  contracting  refers  to  an  amount  or  percent 
by  which  a  contract  price  may  be  adjusted  if  predefined 
contingencies  occur,  such  as  changes  in  the  vendor's  raw 
material  costs  or  ladior  costs.  The  amount  of  the  "escalation" 
is  usually  tied  to  some  predetermined  price  index. 

Synonyms:  Economic  Price  Adjustment 
Antonyms :  None 
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Comments : 

Delete  definition,  add  "A  term  traditionally  used  to  indicate 
the  periodic  price  adjustment  of  a  contract.  It  is  frequently 
computed  by  a  mathematical  formula,  specified  in  the  contract 
or  BOA,  utilizing  well  known  national  indices.  It  is  not 
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uncommon  for  the  contractor  to  be  limited  to  recovery  of  only 
a  portion  of  the  total  fluctuation  defined  by  the  formula  as 
part  of  the  risk  sharing  arrangement  of  the  contract. 

Add  to  (1)  "A  pricing  term". 

Change  (2)  to  "In  Government  contracting  refers  to  an  amount, 
rate  or  percent  by  which  a  contract  price  may  be  adjusted  if 
predefined  contingencies  occur,  such  as  significant  changes 
beyond  its  control  in  the  vendor's  raw  material  costs  or  labor 
costs.  The  amount  of  the  "escalation"  must  be  tied  to  some 
predetermined  price  index." 

Add  to  end  of  first  sentence  in  (1)  "/cost." 

Escalation  would  be  only  an  upward  movement  while  EPA  could  go 
both  ways. 

Eliminate  "more  rarely"  in  (1)  . 

Add  to  (2)  after  "such  as  changes"  -  "upward  or  downward". 

EPA  is  not  synonym  of  Escalation.  EPA  is  application  of 
escalation. 

Should  contract  type  be  included  in  definition,  i.e.  CP  or 
FFP? 

Change  (1)  "application  of  escalation  by  specified  procedures" 
to  "application  of  previously  agreed  price  adjustment  (s)  after 
contract  award. 

Change  (2)  "escalation"  to  "adjustment". 

Change  end  of  (2)  "predetermined  price  index"  to 
"predetermined  public  or  Government  price  index. " 

Change  first  sentence  in  (1)  to  "A  price  increase  or  revision 
upward  due  to  external  influences  such  as  inflation  or  market 
adjustments. " 

Synonyms:  Cost  Growth. 

Antonyms:  Deescalation. 

Revised  Definition: 

ESCALATIOW 

A  pricing  term  traditionally  used  to  indicate  an  upward 
movement  of  price/cost  due  to  inflation  or  market  adjustment. 
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Synonyaui  i  None 
Antonyms :  Deescalation 


Do  you  agree  with  this  definition? 

. 1 . 2 . 3 . . --4 . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 

COMMENT  : 


SYNONYMS : 
ANTONYMS; 
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Original  Definition: 


B3tPBWSB 

Costs  of  operation  and  maintenance  of  activities  on  the 
accrual  basis  for  a  fiscal  period,  as  distinguished  from 
capital  costs  that  will  be  depreciated  over  their  approximate 
service  life. 

Synonyms:  Costs 
Antonyms:  Revenue,  Income 

Survey  Results 


Expense 
First  Round 


Rating  Scale 


Comments : 

Change  definition  to  "A  cost  incurred  in  performance  of  a 
business  operation  or  contract  to  be  accounted  for  on  an 
accrual  basis  for  a  fiscal  year  or  tied  to  a  specific 
contract.  As  distinguished  from  capital  ...." 

Delete  "Revenue,  Income"  from  antonyms 

Change  "will  be  depreciated"  to  "depreciate". 
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Could  break  dovm  definition  into  cash  expenses  and  non-cash 
e3q>enses . 

Change  "on  the  accrual  basis  for  a  fiscal  period,  as  . . . . "  to 
"for  a  fiscal  period." 

Change  definition  to  "The  collection  of  costs  related  to  a 
particular  defined  set  of  activities,  over  a  set  period  of 
time, " 

Change  "Costs  of  operation"  to  "Reasonable  costs,  direct  and 
indirect,  of  operation". 

Why  only  "accrual  basis"?  Is  it  true  for  "cost  basis"? 

Synonyms:  Costs,  Burdens,  Indirect  Costs,  Outgo,  Overhead 

Item,  Consumption,  Spending. 

Antonyms :  Fee . 

Revised  Definition: 

SXPSMSB 

Costs  of  operation  and  maintenance  of  activities  on  the 
accrual  basis  for  a  fiscal  period,  as  distinguished  from 
capital  costs  that  depreciate  over  their  approximate  service 
life. 


Synonyms:  Costs 
Antonymsi  Revenue ,  Income 

Do  you  agree  with  this  definition? 


. 1 . 2--- 

STRONGLY  AGREE 

. 3 . 

AGREE  W/ 

. 4 - 

DISAGREE  W/ 
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Original  Definition! 


Plant  equipment  acquired  by  the  Government,  exceeding  an 
estcd3lished  acquisition  cost  threshold,  used  for  the  purpose 
of  cutting,  abrading,  grinding,  shaping,  forming,  joining, 
testing,  measuring,  heating,  treating  or  otherwise  altering 
the  physical,  electrical  or  chemical  properties  of  materials, 
components  or  other  end  items  entailed  in  manufacturing, 
maintenance,  supply,  processing,  assembly  or  research  and 
development  operations. 

Synonyms :  None 
Antonyms :  None 

Survey  Results 

Industrial  Plant  Equipment 


First  Round 


Rating  Scale 


Comments i 

Add  to  end  of  definition  "anticipated  to  have  value  and  use 
after  the  contract  is  completed." 

Delete  "exceeding  an  established  acquisition  cost  threshold". 
IPE  is  IPE  if  it  is  within  or  in  excess  of  a  cost  threshold. 
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DFARS  245.301. 

Add  to  definition  "This  equipment  is  Government  Furnished 
Equipment  (GFE)  for  the  contractors'  use  in  furtherance  of  the 
Government  contract . " 

Delete  listing  type  of  equipment  -  too  limiting. 

IPE  is  not  necessarily  acquired  by  the  Government  in  all 
cases.  The  contractor  must  sometimes  invest  in  IPE. 

Synonyms :  None 
Antonyms  t  None 

Revised  Definition: 

INDUSTRIAL  PLANT  EQUIPMENT  (IPB) 

Plant  equipment  acquired  by  either  Government  or  industry, 
exceeding  an  established  acquisition  cost  threshold,  used  for 
the  purpose  of  altering  the  physical,  electrical  or  chemical 
properties  of  materials,  components  or  other  end  items 
entailed  in  manufacturing,  maintenance,  supply,  processing, 
assembly  or  research  and  development  operations. 

Synonyms :  None 
Antonyms :  None 

Do  you  agree  with  this  definition? 


. 1 . 2--- 

STRONGLY  AGREE 

. 3 . 

AGREE  W/ 

. 4 . 

DISAGREE  W/ 

- 5 . 

DISAGREE 

. 6 - 

STRONGLY 
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RESERVATION 

DISAGREE 
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SYNONYMS: 

ANTONYMS: 
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Original  Da£lnltlont 

OFgfiBTS 


A  cost  balancing  action  whereby  a  claim  may  be  canceled  or 
lessened  by  a  counterclaim. 

Defective  pricing:  Allowcible  understatements  (e.g., 
counterclaims  or  cost  proposal  errors  that  are  favorable  to 
the  contractor)  which  are  reduced  by  overstatements  of  cost 
that  arise  under  a  defective  pricing  case.  In  order  to 
eliminate  an  Increase  In  the  contract  price  the  offset  cannot 
exceed  the  extent  of  the  overstatement. 

Administrative  Offset:  A  procedure  to  collect  a  debt  owed 
to  the  Government  by  withholding  money  payable  to  contractor 
under  a  contract.  In  order  to  satisfy  the  contractor's  debt 
which  arose  Independently  of  that  contract  and  which  are  in 
compliance  with  the  Federal  Claims  Collection  Act  of  1966. 

Synonyms  t  Counterclaim,  Setoff 
Antonyms  x  None 
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ComMOTtst 

International  offsets  left  out. 

An  offset  may  be  a  deduction  or  credit,  as  well. 

Delete  "In  order  to  eliminate  an  Increase  In  the  contract 
price  the  offset  cannot  exceed  the  extent  of  the 
overstatement . " 

Add  to  paragraph  3  "payable  to  the  contractor  _ " . 

Add  to  paragraph  3  "A  unilateral  procedure  . . . " . 

Add  to  end  of  paragraph  l  "A  tradeoff  wherein  a  cost  Is 
allowed  for  a  particular  segment  of  the  work  but  a 
corresponding  reduction  in  cost  is  agreed  upon  for  another 
segment . " 

Synonyms : 

Antonyms  t 

Revised  Definition: 

QPF8BTS 

A  cost  balancing  action  whereby  a  claim  may  be  canceled  or 
lessened  by  a  counterclaim. 

Defective  pricing:  Allowcd^le  understatements  (e.g., 

counterclaims  or  cost  proposal  errors  that  are  favored^le  to 
the  contractor)  which  are  reduced  by  overstatements  of  cost 
that  arise  under  a  defective  pricing  case.  In  order  to 
eliminate  an  increase  in  the  contract  price  the  offset  cannot 
exceed  the  extent  of  the  overstatement. 

Administrative  Offset:  A  procedure  to  collect  a  debt  owed 
to  the  Government  by  withholding  money  payeUsle  to  the 
contractor  under  a  contract,  in  order  to  satisfy  the 
contractor's  debt  which  arose  independently  of  that  contract 
and  which  are  in  cozr$>liance  with  the  Federal  Claims  Collection 
Act  of  1966. 

Synonyms:  Counterclaim,  Setoff 
Antonyms :  None 


Do  you  agree  with  this  definition? 


. 1 . 2 . 3 . 4 . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 
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CC^IMENT 


SYNONYMS: 

ANTONYMS: 
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Original  Definitions 
WBiqgTBP  AVBWigg  COST  IBTgQP 

A  method  of  determining  the  average  unit  cost  of  inventory  euid 
by  inqplication  an  aid  in  determining  the  cost  of  goods  made, 
sold,  or  held  for  future  sale  or  incorporation  into  higher 
level  end  items.  Under  this  technique,  costs  are  periodically 
con?)uted  by  adding  the  sum  of  the  costs  of  beginning  inventory 
with  the  svun  of  the  costs  of  subsequent  purchases  and  dividing 
by  the  total  niomber  of  units. 

Synonyas  t  None 
Antonyms :  None 
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Survey  Results 

Weighted  Average  Cost  Method 
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Comments t 

Change  "goods  made,  sold,  or  held..."  to  "goods  made  in 
process,  sold  to  regular  customers  or  held...". 

Change  "by  implication  an  aid"  to  "by  historical  conparisons" . 
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Add  to  end  of  definition  "Cost  values  are  obtained  by 
multiplying  the  values  by  their  weights  then  added  together 
and  divided  by  the  sum  of  the  weights . " 

Synonyms:  Unit  Cost  Comparison  Technique. 

Antonyms:  Specific  Identification,  Actual  Cost  Method. 

Revised  Definition: 

WEIGHTED  AVERAGE  COST  METHOD 

A  method  of  determining  the  average  unit  cost  of  inventory  and 
by  implication  an  aid  in  determining  the  cost  of  goods  made, 
sold,  or  held  for  future  sale  or  incorporation  into  higher 
level  end  items.  Under  this  technique,  costs  are  periodically 
confuted  by  adding  the  sum  of  the  costs  of  beginning  inventory 
with  the  siun  of  the  costs  of  subsequent  purchases  and  dividing 
by  the  total  number  of  units. 

Synonyms:  Unit  Cost  Comparison  Tscbnlque 
Antonyms :  None 

Do  you  agree  with  this  definition? 

. 1 . 2 . 3 . 4 . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 

C0I®1ENT: _ 


SYNONYMS: 

ANTONYMS: 
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Original  Dafinitiont 
AllQgfttiQB 


Funding;  An  amount  of  money  treuisf erred  from  one  agency, 
bureau  or  account  that  is  set  aside  in  an  appropriation  of  the 
various  committees  having  spending  responsibilities  to  carry 
out  the  purposes  of  the  parent  appropriation  or  fund. 

Within  DOD,  the  money  is  being  transferred  from  the  services 
to  the  appropriate  MAJCOMS. 

Financial;  A  cost  accounting  procedure  which  results  in  a 
reasoncdile  distribution  of  costs  among  one  or  more  cost 
objectives  (e.g.,  products,  programs,  contracts,  and 
activities)  .  This  includes  both  direct  assignment  of  costs  and 
the  reassignment  of  a  share  from  an  indirect  pool. 

Synonyms :  None 
Antonyms :  None 

Survey  Results 

A I  location 


First  Round 


Rating  Scale 


L3 


Comnents  t 


Perhaps  change  reasonable  distribution  to  appropriate 
distribution. 

Funding:  An  amount  of  money  distributed  or  assigned  by  formal 
action  to  a  particular  group  or  account  for  a  particular  use 
or  period  of  time. 

Addition  of  the  meaning  of  term  for  resources  in  the 
acquisition  arena,  i.e.,  allocation  of  manpower,  skills, 
computer  time  and  memory  to  tasks. 

Add  "Represents  money  that  can  be  obligated."  to  Funding. 

Delete  from  Funding  -  "transferred  from  one  agency,  bureau  or 
account " . 

Delete  "Within  DOD"  paragraph. 

Change  "from"  to  "to"  and  "parent  appropriation  or  fund."  to 
"concerned  organization"  in  Funding  paragraph. 

Add  "from  DOD  to  the  services"  to  "Within  DOD"  paragraph. 

Change  "transferred"  to  "earmarked  for"  in  Funding  paragraph. 

Change  "objectives"  to  "categories"  in  Financial  paragraph. 
We  are  not  setting  objectives  when  we  allocate  costs. 

Change  "one  agency,  bureau  or  account"  to  "entity"  in  Funding 
paragraph . 

Not  sure  that  funds  had  to  be  transferred  -  allocation  could 
be  done  by  notation  or  journal  entry. 

Delete  "of  the  various  committees  having  spending 
responsibilities"  from  Funding  paragraph. 

Change  "transferred  from  the  services  to  the  appropriate 
MAJCOMS."  to  "flowed  down  from  higher  headquarters  to  the 
appropriate  users . " 

Change  Financial  to  read  "An  accounting  procedure  which 
assigns  costs  to  an  identified  usage  or  purpose." 

Change  "agency,  bureau  or  account"  to  "e.g.  agency  bureau  or 
account"  to  not  limit  definition. 

Spell  out  acronym. 
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Change  Financial  to  read  "A  cost  accounting  process  of 
assigning  a  cost,  or  group  of  costs,  to  one  or  more  cost 
objectives,  in  reasonable  and  realistic  proportion  to  the 
benefit  provided  or  other  equitable  relationship. 

Delete  "Within  DOD"  paragraph. 

Consider  the  Accounting  definition  "A  systematic  distribution 
or  assignment  of  a  total  amount  among  several  years,  accounts, 
products ,  departments  or  other  elements . " 

Synonyms:  Allotment,  earmark,  assignment,  allowance,  portion, 
quota,  share  allotment,  set  aside. 

Antonyms:  Double  Counting. 


Revised  Definition: 

Allocation 

Funding;  An  amount  of  money  In  a  Government  appropriation 
transferred  to  an  agency,  bureau  or  account  having  spending 
authority  to  carry  out  the  purposes  of  the  that  appropriation. 


Financial;  A  cost  accounting  procedure  which  results  in  a 
reasonable  distribution  of  costs  among  one  or  more  cost 
objectives  Including  products,  programs,  contracts,  and 
activities.  This  includes  both  direct  assignment  of  costs  and 
the  reassignment  of  a  share  from  an  indirect  pool. 

Synonyms :  None 
Antonyms :  None 

Do  you  agree  with  this  definition? 

. 1 . 2 . 3 . ---4 . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 

COMMENT; 


SYNONYMS: 

ANTONYMS: 
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Original  Dafinition: 


The  systematic  reduction  of  an  indebtedness  or  recorded  asset 
value  over  a  specific  period  of  time  by  periodic  payments  to 
a  creditor  or  charges  to  an  expense,  in  accordance  with 
generally  accepted  accounting  procedures  or  principles. 

Synonyms t  Liquidation,  Allocation,  Writeoff 
Antonyms*  Direct  charge 

Survey  Results 
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Comments : 

Reference  to  GAAP  not  necessary. 

Delete  allocation  as  synonym. 

Change  "asset  value"  to  "asset  net  value". 
Delete  direct  charge  as  an  antonym. 

Change  "reduction"  to  "liquidation". 
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Change  "charges  to  an  expense”  to  "charges  against  a  capital 
account " . 

Change  definition  to  "A  system  or  method  which  reflects  how 
much  of  the  value  of  an  asset  is  reduced  due  to  usage  or  the 
passage  of  time." 

Add  "Amortization  is  often  calculated  to  occur  over  a 
specified  period  of  time." 

Add  after  "asset  value",  "usually  a  depreciable  capital 
asset" . 

Delete  "recorded"  and  "value"  from  "recorded  asset  value". 
Add  after  "expense",  "account". 

Change  "reduction  of  an  indebtedness"  to  "extinguishment  of 
debt " . 

Synonyms:  Depreciation  (?) ,  reduction,  reden5)tion. 

Antonyms:  Expense  item.  Expensed. 

Revised  Definition: 

Amortization 

The  systematic  liquidation  of  an  indebtedness  or  recorded 
asset  value  over  a  specific  period  of  time  by  periodic 
payments  to  a  creditor  or  charges  to  an  expense  account,  in 
accordance  with  generally  accepted  accounting  procedures  or 
principles . 

Synonyms:  Liquidation,  Writeoff 

Antonyms :  None 

Do  you  agree  with  this  definition? 

. 1 . 2 . 3 . . 4 . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 

COMMENT; 


SYNONYMS  ; 
ANTONYMS; 
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Original  Definition t 
Bid 

An  offer  to  perform  a  contract  by  providing  lessor  and  or 
material  for  a  specific  price.  In  federal  government 
contracting,  this  offer  is  provided  in  response  to  an 
invitation  for  bid. 


Synonyms :  None 
Antonyms  t  None 

Survey  Results 


U) 

0) 

Ul 

c 

o 

a 

ui 

« 

q: 


v 

c 

(U 

u 

L 

0) 

0. 


Rating  Scale 


Comments  x 

Change  "specific  price"  to  "fixed  or  specific  price". 

Axe  "labor  and  or  material"  too  specific  and  do  they  include 
special  test  equipment/technical  services. 


Change  "perform  a  contract  by  providing  IcUoor  and  or  material" 
to  "provide  supplies  or  services  for  a  contract". 


Capitalize  "Invitation  for  Bid". 


Change  "Bid"  to  "Proposal"  -  outdated  term. 

Add  to  first  sentence,  "in  accordance  with  specified  contract 
terms. " . 

Reference  legal  obligation  in  definition. 

Change  "and  or  material"  to  "and/or  material". 

Change  "perform"  to  "form". 

Change  "invitation  for  bid"  to  "a  solicitation  which  invites 
the  submission  of  such  offers.". 

Address  "sealed  bid". 

Add  to  first  sentence,  "A  firm- fixed-price,  usually 

irrevoc2d)le,  offer  to  perform  ...". 

Add  to  end  of  definition,  "that  will  not  be  negotiated.". 

Add  to  end  of  first  sentence,  "usually  on  a  firm  fixed  price 
basis. " 

Change  second  sentence  to  "To  perform  the  work  specified  in  an 
invitation  for  bid  (IFB)." 

Change  first  sentence  to  "An  offer  to  perform  the  scope  of 
work  specified  in  a  contract  for  a  specific  price." 

Don't  limit  definition  to  Government  only. 

Change  first  sentence  to  "An  offer  by  a  prospective  purchaser 
to  buy  goods  or  services  at  a  stated  price,  or  an  offer  by  a 
prospective  seller  to  sell  his  goods  or  services  for  a  stated 
price. 

Synonyms  t  Offer . 

Antonyms:  Request  for  Proposals  (RFP) 

Revised  Definition: 

Bid 

An  offer  to  perform  a  contract  by  providing  goods  or  services 
for  a  specific  price.  In  Federal  Government  contracting,  it 

is  the  technical  term  for  an  Irrevocable  offer  in  response  to 
an  Invitation  For  Bid  (IFB) . 
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Synonyms t  Offer,  Proposal. 
Antonyms  t  None 


Do  you  agree  with  this  definition? 


. 1 . 2 . 3 - 4 . 5 . . -6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 


COMMENT: 


SYNONYMS: 

ANTONYMS: 
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Origloal  Deflnltlont 

The  act  by  an  authorized  individual  affirming  the  intent  of  an 
agency  or  company  to  take  or  accept  a  defined  action  not  yet 
formalized  by  execution  of  a  contract. 

Funding:  A  firm  administrative  reservation  of  funds  based  upon 
firm  procurement  directions,  orders,  requisitions,  certified 
purchase  requests,  and  budgetary  authorizations  which  set 
aside  certain  funds  for  a  particular  contract  without  further 
recourse  to  the  official  responsible  for  certifying  the 
availability  of  funds. 

Within  DOD,  reservation  of  funds  are  set  aside  by  the 
appropriate  operating  division  (wing  or  base)  for  use  on  a 
particular  item. 

Accounting:  The  method  of  accounting  for  the  availcdsle  balance 
of  an  appropriation,  fund,  or  contract  authorization  whereby 
commitments  are  recorded  in  the  accounts  as  reductions  of  the 
availedale  balance. 

Synonyms  t  None 
Antonyms :  None 

Survey  Results 


Comm  I tment 

First  Round 


A  im  0  to 

Rating  Scale 
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Coomants  t 


Change  first  paragraph  to  "An  engagement  by  contract  involving 
financial  obligation." 

Change  last  paragraph  to  "Accounting:  The  total  accumulated 
financial  obligation  against  a  contract  or  other  objective  at 
a  specific  point  in  time." 

Note:  A  commitment  can  be  either  written  or  unwritten. 

Second  paragraph  add  "It  indicates  intention (s)  to  incur 
obligations. " 

Third  paragraph  change  "use  on  a  particular  item. "  to 
"something  to  be  bought  in  the  future." 

Third  paragraph  delete  " (wing  or  base) " . 

Conflict  between  first  and  second  paragraph  -  "affirming  the 
intent"  and  "A  firm  reservation  of  funds  based  on  firm 
procurement  directions". 

Third  paragraph  change  "division  (wing  or  base) "  to 
"organizations" . 

Second  paragraph  change  "contract"  to  "activity".  Funding  can 
be  set  aside  for  other  than  contracts,  i.e.  interagency 
agreements . 

Delete  first  paragraph  -  intent  is  not  binding. 

First  paragraph  change  "agency  or  company"  to  "entity". 
Paragraph  three  delete  "reservation  of". 

Second  paragraph  change  "firm  administrative"  to  "definitive. 

Second  paragraph  change  "firm  procurement  directions"  to 
"procurement  directives". 

Synonyms :  None 
Antonyms  t  None 

Revised  Definition! 

Commitment 

The  act  by  an  authorized  individual  affirming  the  intent  of  an 
agency  or  company  to  take  or  accept  a  defined  action  not  yet 
formalized  by  execution  of  a  contract. 
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Funding:  An  administrative  reservation  of  funds  based  upon 
procurement  directions,  orders,  requisitions,  certified 
purchase  requests,  and  budgetary  authorizations  which  set 
aside  certain  funds  for  a  particular  contract  without  further 
recourse  to  the  official  responsible  for  certifying  the 
availability  of  funds. 

Within  DOD,  funds  are  set  aside  by  the  appropriate  operating 
organizations  for  use  on  a  particular  item. 

Accounting:  The  method  of  accounting  for  the  availcdile  balance 
of  an  appropriation,  fund,  or  contract  authorization  whereby 
commitments  are  recorded  in  the  accounts  as  reductions  of  the 
available  balance. 

Synonyms :  None 
Antonyms  t  None 

Do  you  agree  with  this  definition? 

. 1 . 2 . 3- . 4 .  -5 .  -6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 


COMMENT: 


SYNONYMS: 


ANTONYMS: 
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Original  Definition: 
Cost 


(1)  For  the  Seller:  The  amount  of  money  or  equivalent 

incurred  for  supplies  or  services  exclusive  of  profit  or  fee. 

(2)  For  the  Buyer;  The  amount  of  money  or  equivalent  paid 
for  supplies  or  services  including  profit  or  fee. 

Synonyms:  Expense,  Consideration,  Charge,  Total  Cost 
Antonyms :  None 


Survey  Results 
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Comments : 

Change  (l)  to  "The  total  eunount  of  money  or  equivalent 
incurred  for  the  production  or  purchase  of  supplies  or  the 
performance  of  services  exclusive  of  profit  or  fee." 

Change  (2)  to  "The  amount  of  money  or  equivalent  paid  for 
supplies  or  services  including  the  seller's  profit  or  fee,  the 
seller's  price." 


Remove  "Total  Cost"  from  Synonyms  because  total  cost  may 
Include  fee. 


Remove  "Consideration"  from  Synonyms. 

Add  (3)  "A  direct  or  indirect  charge  of  a  specific  or  unique 
element  allocated  to  a  particular  cost  objective." 

Consider  type  of  contract  -  Cost  Plus  Fixed  Fee,  cost  is 
separate  from  profit. 

Synonyms:  Direct  Cost,  Indirect  Cost,  Billed  Amount,  Actual 
Cost. 

Antonyms:  Applicable  Credit,  Negative  E3q>enditure . 

Revised  Definitions: 

Cost 

(1)  For  the  Seller:  The  amount  of  money  or  equivalent 

incurred  for  supplies  or  sezrvices  exclusive  of  profit  or  fee. 

(2)  For  the  Buyer;  The  amount  of  money  or  equivalent  paid 
for  supplies  or  services  including  the  seller's  profit  or  fee. 

Synonyms:  Expense,  Consideration,  Charge,  Total  Cost 
Antonyms :  None 

Do  you  agree  with  this  definition? 

. 1 . 2 . 3 . 4 . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  H/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 

COMMENT;  _ 


SYNONYMS: 

ANTONYMS: 
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Original  Definition: 


Dellnonencv 

(1)  Failure,  omission,  or  violation  of  contractual 
obligation  or  duty. 

(2)  The  actual  failure  by  the  contractor  to  meet  the 
contract  delivery  or  performance  schedule,  or  the  potential 
failure  to  do  so  by  failing  to  maintain  required  progress 
in  contract  performance  as  required  by  the  contract  delivery 
or  performance  schedule 


Synonyms:  Overdue,  Tardy,  Late 
Antonyms:  Early,  Accelerated,  Timely 

Survey  Results 
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Rating  Scale 


Comments: 

In  (2)  change  "required  progress"  to  "progress". 

Change  (2)  to  "The  actual  failure  by  a  contractor, 
subcontractor  or  supplier  to  meet  the  contract  delivery  or 
performance  schedule,  or  the  potential  failure  to  do  so  by  not 


maintaining  progress  as  required  by  the  contract  delivery  or 
performance  schedule." 


Delete  (2) . 

Is  "potential  delinquency"  a  delinquency? 

In  (2)  delete  "potential". 

In  (2)  add  Government  contribution  to  delinquency  by  failing 
to  deliver  GFE  on  time. 

Change  (2)  from  "meet  the  contract  delivery  or  performance 
schedule"  to  "meet  the  contract  delivery  or  performance 
schedule  or  performance  requirements". 

Synonyms:  Pass  Due,  Deficient,  Substandard  Performance, 

Breach,  Noncompliance,  In  Default,  Derelict,  Failure, 
Behind  Schedule,  Missed  Milestone. 

Antonyms:  Proficient,  Standard  Performance,  Compliant. 

Revised  Definition: 

PtliagW«tfiY 

(1)  Failure,  omission,  or  violation  of  contractual 

obligation  or  duty. 

(2)  The  actual  failure  by  the  contractor  to  meet  the 
contract  delivery  or  performance  schedule,  performance 
requirements  or  by  failing  to  maintain  required  progress  in 
contract  performance  as  required  by  the  contract  delivery  or 
performance  schedule 

Synonyms:  Overdue,  Tardy,  Late 
Antonyms:  Early,  Accelerated,  Timely 

Do  you  agree  with  this  definition? 

. 1 . 2 . 3 . 4 . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 

AGREE  RESERVATION  RESERVATION  DISAGREE 

COMMENT:  _ 


SYNONYMS: 

ANTONYMS: 
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Original  Daflnltlont 


Executed  Contract 

A  written  dociunent  which  has  been  signed  by  both  parties  and 
mailed  or  otherwise  furnished  to  each  party,  which  expresses 
the  requirements,  terms,  and  conditions  to  be  met  by  each 
party . 


SynonyiBS :  None 
Antonyms :  None 

Survey  Results 


Executed  Contract 

First  Round 

XX  I - 
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Rating  Scale 


Comments : 

Change  "A  written  document  which. . . "  to  "A  written  docximent  in 
the  hands  and  under  control  of  all  parties  which..." 

Executed  or  Executory  Contract. 

Change  "and  mailed  or  otherwise"  to  "and". 

Should  price/consideration  be  added? 


Change  "which  expresses"  to  "which  clearly  expresses  the 
mutually  agreed" . 

Does  definition  cover  new  technology,  i.e.  EDI,  FAX,  etc.? 


Synonymsi  Covenant,  Legally  Binding  Agreement,  Definitized 
Contract,  Award,  Purchase  Agreement,  Fully  Signed 
Document . 

Antonyms*  Ratified,  Executory  Contract. 

Revised  Definition: 

Executed  Contract 

A  written  document  which  has  been  signed  by  both  parties  and 
furnished  to  each  party,  which  expresses  the  requirements, 
terms,  and  conditions  to  be  met  by  each  party. 


Synonyms:  Definitized  Contract 
Antonyms :  None 

Do  you  agree  with  this  definition? 

. 1 . 2 . 3 . 4 . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 


COMMENT: 


SYNONYMS: 

ANTONYMS: 
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Original  Definition: 

Independent  Co*«- 

A  cost  estimate  developed  outside  the  normal  advocacy 
channels,  independent  of  any  cost  information  provided  by  the 
offeror,  used  for  the  purpose  of  comparing  with  bids  or 
proposals.  Preparation  of  independent  costs  estimates 
generally  include  representations  from  the  areas  of  cost 
analysis,  procurement,  production  management,  engineering,  and 
program  management. 

Synonyms:  Independent  Government  Cost  Estimate  (IGCE) 
Antonyms :  None 


Survey  Results 
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Rating  Scale 


Comments : 

Change  "cost  analysis"  to  "cost  estimating". 

Change  "include  representations"  to  "includes  input". 
Eled>orate  on  advocacy  channels. 
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"within  the  procuring  organization" . 


Add  after  "developed"  - 


Add  after  "proposals"  -  "and  often  used  in  negotiations." 

Add  after  "representations  from"  -  "one  or  more  areas  of 
price/cost  analysis". 

Put  parenthesis  around  sentence  2. 

Synonymst  Should  Cost  Estimate. 

Antonyms I  Contractor  Prepared  Cost  Information,  Dependent 
Cost  Estimate. 

Revised  Definition t 

Independent  Coat  Batiaatc 

A  cost  estimate  developed  independent  of  any  cost  information 
provided  by  the  offeror,  used  for  the  purpose  of  comparing 
with  bids  or  proposals.  Preparation  of  independent  costs 
estimates  generally  includes  representations  from  one  or  more 
of  the  areas  of  cost/price  analysis,  procurement,  production 
management,  engineering,  and  prograun  management. 

Synonyms:  Independent  Government  Cost  Estimate  (IGCE) 

Antonyms:  Contractor  Prepared  Cost  Information 

Do  you  agree  with  this  definition? 

. 1 . 2 . 3 . 4 . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 

CC»®4ENT: 


SYNONYMS: 


ANTONYMS: 
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Original  Dafinition: 


Novation  Aoraamant 

A  legal  instrmnent,  executed  by  the  parties  to  a  contract  and 
a  successor  in  interest,  which  transfers  all  obligations  and 
rights  under  the  contract  to  the  successor. 

The  government  may  recognize  a  third  party  as  a  successor  of 
a  government  contract  when  the  third  party's  interests  arises 
out  of  the  transfer  of  1)  all  the  contractor's  assets,  or  2) 
the  entire  portion  of  the  assets  involved  in  the  performing  a 
contract . 

Synonyms :  None 
Antonyms  t  None 

Survey  Results 
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Comments : 

Change  "performing"  to  "performance  of". 

Change  "under  the  contract  to  the  successor"  to  "under  the 
contract  of  one  party  to  its  successor  in  interest . " 
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Paragraph  2  change  "successor"  to  "successor  in  interest". 

Novation  agreements  can  be  made  only  when  a  con^jany  changes 
their  name  only.  (?) 

Paragraph  2  change  "may"  to  "reserves  the  right  to". 

Synonyms:  Transfer  Agreement,  Discharge  of  Contract,  Mutual 
Rescission,  Cancellation,  Substituted  Contract, 
Contract  Name  Change. 

Antonyms : 

Revised  Definition: 

Novation  Agreement 

A  legal  instrument,  executed  by  the  parties  to  a  contract  and 
a  successor  in  interest,  which  transfers  all  obligations  and 
rights  under  the  contract  of  one  party  to  the  successor  In 

Interest. 

The  government  reserves  the  right  to  recognize  or  not 
recognize  a  third  party  as  a  successor  In  Interest  of  a 
government  contract  when  the  third  party's  interests  arises 
out  of  the  transfer  of  1)  all  the  contractor's  assets,  or  2) 
the  entire  portion  of  the  assets  involved  in  the  performance 
of  a  contract. 

Synonyms :  None 
Antonyms:  None 

Do  you  agree  with  this  definition? 

- 1 . 2 . ---3- . ...4 . 5-. . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 

COMMENT: 


SY'PIONYMS : 
ANTONYMS: 
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Original  Definition: 
RQYftltY 


Condensation  paid  to  the  owner,  vendor  or  lessor  of  personal, 
real,  tangible  or  intangible  property  for  the  use  of  that 
property.  Usually  a  percentage  of  the  selling  price  of  goods 
and  services,  production  of  which  employs  the  property. 

Synonyms:  Commission  Payment,  Use  Fee 
Antonyms:  Royalty  Free  Use 

Survey  Results 

Roya I ty 


First  Round 


SA  A  AR  DR  D  so 


Rating  Scale 


Comments : 

Delete  "production  of  which  endloys  the  property" . 

Sentence  2  not  clear. 

Do  not  consider  a  fee  for  rental  of  property  to  be  royalty. 
Protected  rights  of  the  owner. 

Synonyms:  Intellectual  Property  Fee,  License  Fee,  Rent. 

Antonyms:  No  Fee,  Rent  Free. 
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Revised  Definition} 


Rgyftlty 

Compensation  paid  to  the  owner  or  vendor  of  personal,  real, 
tangible  or  intangible  property  for  the  use  of  that  property. 


Synonyms t  Commission  Payment,  Use  Fee,  Intellectual  Property 
Fee. 

Antonyms:  Royalty  Free  Use 

Do  you  agree  with  this  definition? 

. 1 . 2 . 3 - 4 . 5 . 6 - 

STRONGLY  AGREE  AGREE  W/  DISAGREE  W/  DISAGREE  STRONGLY 
AGREE  RESERVATION  RESERVATION  DISAGREE 


COMMENT: 


SYNONYMS; 

ANTONYMS; 
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